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I.  OVERVIEW 

A.  SUMMARY 

The  Arizona  Center  for  Mathematical  Sciences  (ACMS)  was  formed  in  1986  with  an  initial  three- year  grant 
from  the  Air  Force  Office  of  Scientific  Research  under  the  University  Research  Initiative  Program.  Offices 
and  computer  facilities  for  ACMS  personnel,  originally  located  on  the  eighth  floor  of  the  Gould-Simpson 
Building,  were  moved  last  summer  to  a  larger  space  in  the  Economics  Building,  Close  ties  are  maintained 
with  faculty  in  the  Departments  of  Mathematics,  Physics,  Aerospace  and  Mechanical  Engineering,  Optical 
Sciences  and  the  Applied  Mathematics  Program  at  the  University  of  Arizona. 

The  ongoing  and  primary  goal  of  current  research  is  the  pursuit  of  understanding  of  nonlinear  processes 
in  natural  phenomena  arising  in  optics.  Optics,  a  fertile  area  of  research  in  nonlinear  science,  provides 
a  vital  underpinning  for  new  technology  and  also  serves  as  a  useful  paradigm  for  gaining  an  increased 
understanding  in  other  fields.  For  example,  turbulence  in  optics,  the  study  of  the  complex  space-time 
filaments,  patterns  and  defects  which  appear  in  wide  aperture  high-power  lasers  and  counterpropagating 
beams,  may  be  more  analytically  tractable  than  in  other  branches  of  continuum  mechanics. 

This  report  details  research  progress  made  under  AFOSR  contract  AFOSR-F49620-94-1-0144DEF,  “Semi¬ 
conductor  Laser  Dynamics”  during  the  period  of  the  contract.  As  specified  in  the  original  proposal,  our 
overall  goal  was  to  establish  a  close  synergism  between  sophisticated  mathematical  analysis,  fundamental 
physics  and  high  performance  computing  with  interactive  graphical  visualization.  An  interdisciplinary 
team  of  researchers  from  mathematics,  optical  sciences  and  physics  has  been  assembled  locally  at  the  Uni¬ 
versity  of  Arizona,  close  collaborative  ties  are  being  maintained  with  Air  Force  Scientists  at  the  Philfips 
Laboratory,  Kirtland  AFB  and  with  experts  at  outside  institutions.  We  have  recently  established  a  di¬ 
rect  working  link  with  Opto  Power  Corporation,  a  semiconductor  laser  company  that  moved  to  Tucson 
in  1996.  Dissemination  of  research  results  is  facilitated  through  regular  ACMS  Workshops,  national  and 
international  conference  invited  and  contributed  papers,  regular  laboratory  visits,  email,  fax  and  telephone 
communications . 

B.  COLLABORATIONS 

A  unique  strength  of  ACMS  activity,  is  the  degree  to  which  collaboration  is  fostered  between  local  re¬ 
searchers  at  the  University  of  Arizona,  with  research  scientists  at  the  Air  Force  laboratories  and  through 
international  links  to  the  European  Community.  Regular  visits  are  made  to  the  Phillips  Laboratory,  Kirt¬ 
land  AFB,  Albuquerque  for  the  express  purpose  of  disseminating  research  results  and  obtaining  direct 
feedback  from  Air  Force  scientists  on  our  work.  Our  direct  contacts  at  the  Laboratory,  Dr.  D.  Depatie, 
Drs.  D.  Bossert  and  M.M,  Wright,  are  in  regular  contact  by  phone,  fax  and  email.  One  of  the  ACMS 
graduate  students,  Mr.  J.K.  White  spent  a  summer  period  workiing  closely  with  D,  Bossert  on  the  study 
of  filamentation  instabilities  in  broad  area  high-power  semiconductor  lasers.  A  number  of  joint  papers  with 
this  group  have  been  published  or  recently  submitted. 

The  Center  has  acted  as  the  host  for  annual  AFOSR  contract  review  workshops  which  bring  together 
groups  funded  directly  by  AFOSR  to  meet  and  interact  with  Air  Force  scientists.  The  last  such  meeting 
“AFOSR/ACMS  Nonlinear  Optics  Workshop”  was  held  from  October  10-12.  1996. 

A  direct  industrial  link  with  Opto  Power  Corporation,  a  manufacturer  of  high  power  semiconductor  laser 
sources  was  established  in  May  1996.  We  have  been  working  with  this  company  to  design  high  brightness 
sources  using  our  laser  simulation  capability. 

Our  activity  on  the  international  front  offers  a  unique  window  on  relevant  European  science  and  technology. 
Two  years  ago,  an  international  workshop  on  “Singularities  in  Patterns  and  Collapse:  Applications  to 
Semiconductor  Lasers  and  Critical  Focusing  of  Ultra-Short  Optical  Pulses”,  was  co-organized  by  the  present 
PI  and  by  Professor  J.G.  Mclnerney,  Professor  of  Physics  at  University  College  Cork,  Ireland  and  a 
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collaborator  on  the  present  contract.  This  workshop,  sponsored  mainly  by  BOARD  (London),  was  a 
resounding  success,  drawing  together  Air  Force  scientists  and  leading  academic  and  industrial  experts 
from  Europe  and  the  U.S.  The  leading  European  scientists  attending  the  meeting  were  participants  in 
two  European  Union  semiconductor  laser  network  contracts.  A  follow-on  BOARD- sponsored  Workshop 
entitled  “Fundamentals  and  Modeling  of  Lasers  and  Ultrashort  Pulse  Interactions”,  will  be  held  from  July 
20  -  26,  1997. 

C.  RESEARCH  PROGRAM 

The  overall  research  program  is  succinctly  summarized  in  the  block  diagrams  of  Figs.  1  and  2.  These  list 
the  research  methodology  (Fig.  1)  and  individual  research  projects  (Fig.  2).  Developing  robust  mathemat¬ 
ical  and  physical  models  to  describe  semiconductor  amplifiers  and  lasers,  requires  detailed  knowledge  of  the 
microscopic  many-body  physics  of  the  interaction  of  light  with  a  semiconductor  material.  At  the  same  time 
one  needs  an  understanding  of  the  methodology  used  in  the  reduction  of  realistic  complex  mathematical 
models  of  a  physical  problem  to  more  manageable  universal  order  parameter  equation  descriptions.  Thus 
our  pooled  expertise  simultaneously  brings  to  bear  fundamental  solid  state  theoretical  physics  with  asymp¬ 
totic  and  singular  perturbation  methods  in  mathematics.  Computation  links  these  two  diverse  disciplines 
and  is  carried  out  on  aU  levels,  ranging  from  pre-  and  post-  processing  on  high  performance  graphics  work¬ 
stations  through  to  intermediate  size  in-house  supercomputing  and  on  to  remote  supercomputing  on  vector 
and  massively  parallel  machines  (CM-5  at  AHRPC  (Minnesota),  SP-2  (Maui)  and  Cray  C-90  (WES).  Our 
recent  acquisition  of  a  Silicon  Graphics  Power  Challenge  L  4-processor  compute  engine,  an  Indigo-2  and  4 
Indy  graphics  workstations  under  the  recent  DURIP  program  has  greatly  enhanced  our  in-house  computing 
capability. 

The  diversity  of  our  overall  research  program  reflects  the  committment  made  in  the  original  research 
proposal  to  unify  fundamental  physics  and  mathematics  in  order  to  establish  a  hierarchy  of  robust  math¬ 
ematical  models  capable  of  addressing  the  important  application  area  of  high  power  semiconductor  laser 
sources.  The  following  itemized  summaries  offer  a  global  picture  of  our  research  program  and  indicate  how 
the  individual  parts  contribute  to  achieving  the  specific  long  term  goal. 

•  Femtosecond  Pulse  Propagation  in  Semiconductor  Amplifiers  and  Lasers 

The  study  of  femtosecond  pulse  propagation  in  semiconductor  amplifiers  (narrow  stripe  and  broad 
area)  is  allowing  us  to  establish  the  fundamental  physics  responsible  for  limiting  the  actual  modula¬ 
tion  bandwidth  achievable  in  semiconductor  lasers  or  amplifiers.  Dynamic  bandgap  renormalization, 
spectral  hole  burning,  Coulomb  screening,  plasma  cooling  and  non-MarkolRan  dynamics  aU  signifi¬ 
cantly  influence  the  details  of  pulse  propagation  in  semiconductor  amplifiers. 

•  MicroscopicaRy  Computed  Quantum  Well  Gain  Spectra 

An  important  open  question  from  the  outset,  has  been  how  to  effect  a  closure  of  the  infinite  hierarchy 
of  many-body  equations  describing  the  semiconductor  optical  response.  A  significant  breakthrough 
was  made  by  us  recently  in  resolving  the  significant  discrepancy  between  experimentaRy  measured 
and  theoreticaRy  computed  Quantum  WeR  gain  spectra.  Although  exceRent  agreement  had  been 
achieved  with  bulk  gain  spectra,  it  was  found  necessary  to  include  the  fuR  nondiagonal  scattering 
terms  in  the  Quantum  Boltzmann  equation  to  achieve  quantitative  agreement  with  Quantum  Well 
structures. 

•  Universal  Order  Parameter  Equation  Description  of  Laser  Patterns 

The  ambitious  goal  of  deriving  a  robust  rate  equation  description  of  wide  aperture  semiconductor 
lasers  offers  a  huge  payoff  in  terms  of  modeling  of  real  structures.  The  current  state  of  the  art 


2 


modeling  of  wide  aperture  lasers  and  amplifiers  uses  the  Beam  Propagation  Method  (BPM)  and  does 
not  account  for  time-dependent  behavior.  On  the  other  hand,  the  full  Maxwell-Semiconductor  Bloch 
many-body  laser  equcitions  are  weU  beyond  the  capabilities  of  the  most  powerful  massively  parallel 
or  vector  supercomputers  available  today  or  in  the  foreseeable  future.  By  working  with  simpler 
prototypical  laser  .models,  such  as  the  two-level  laser  MaxweU-Bloch  equations,  we  can  evaluate  the 
robustness  of  the  order  parameter  equation  description  against  the  original  physical  model.  This 
experience  is  essential  as  a  precursor  to  tackling  the  fuU  semiconductor  laser  equations.  We  have 
recently  derived  a  Semiconductor  Laser  Complex  Swift-Hohenberg  equation  description  which  has 
been  verified  against  a  full  scale  nonlinear  pde  model. 

•  Controlling  Optical  Turbulence 

Our  work  under  the  project  “Universal  Order  Parameter  Equation  Description  of  Laser  Patterns” 
provided  the  impetus  for  extending  the  complex  order  parameter  equation  description  to  both  single 
longitudinal  and  two  longitudinal  mode  wide  aperture  semiconductor  lasers.  The  “finite  support 
quasi-plane  traveling  wave”  solutions  of  the  Complex  Swift  Hohenberg  equations  were  shown  to  be 
ubiquitous  for  all  classes  of  wide  aperture  lasers.  This  motivated  us  to  implement  a  feedback  control 
scheme  which  simultaneously  filters  spatially  and  spectrally  the  weakly  turbulent  output  of  a  broad 
area  laser.  The  scheme  has  been  demonstrated  to  be  extremely  robust,  stabilizing  and  steering  the 
laser  output. 

•  Modeling  Spatio-Temporal  Dynamics  of  Broad  Area/Flared  Semiconductor  Laser  Amplifiers 

This  project  has  a  direct  bearing  on  the  applications  areas  of  direct  interest  to  the  Phillips  Laboratory 
scientists  and  the  commercial  needs  of  Opto  Power  Corporation.  A  rate  equation  approach  incorpo¬ 
rating  the  microscopic  many-body  physics  through  a  look-up  table  for  the  gain  and  refractive  index 
functions,  allows  in  principle,  for  a  full  study  of  the  spatiotemporal  evolution  of  filaments  in  broad 
area/flared  amplifiers  and  lasers.  Simulations  of  laser  switch-on  dynamics,  filamentation  instabilities 
in  broad  area  lasers  and  flared  MOPA  structures,  feedback  instabilities  of  MOPA’s  as  a  function  of 
feedback  level  and  eventually,  nonlinear  stabilization  and  control  methods  for  such  structures  are  all 
within  the  capabilities  of  this  level  of  modeling.  Problems  with  spurious  high-transverse  wavenumber 
instabilities  were  taken  care  of  by  adding  an  artificial  diffusion  term  to  the  nonlinear  pde’s  describing 
the  laser.  A  significant  improvement  on  this  model  is  discussed  next. 

•  Full  Space-Time  Simulation  of  High  Brightness  Semiconductor  Lasers 

Significant  progress  has  been  made  in  constructing  a  mathematical  model  of  wide  aperture  semi¬ 
conductor  lasers  which  simultaneously  resolves  multi-longitudinal  mode  and  transverse  filamentation 
instabilities  in  high  brightness  lasers.  The  original  working  model  discussed  under  project  “Simula¬ 
tions  of  Broad  Area  Flared  Amplifiers  and  Lasers”  utilized  a  nonlinear  gain/index  model  but  ignored 
gain  and  index  dispersion.  Gain  and  index  spectra  computed  microscopically  for  specific  laser  struc¬ 
tures  and  including  the  contributions  of  barrier  states,  are  fed  as  rational  functions  approximations 
into  the  laser  simulation  code.  Using  this  model,  we  have  been  able,  for  the  first  time,  to  simulate  the 
full  space-time  evolution  of  all  the  relevant  optical  fields  and  carrier  density  in  the  state-of-the-art 
commercial  MOPA  high  brightness  source.  Our  simulations,  suggesting  a  redesign  of  this  device  for 
improved  stability,  has  attracted  the  attention  of  Opto  Power  Corporation  laser  design  engineers. 

•  Plasma  and  Lattice  Heating  Effects  in  Quantum  Well  Lasers  and  VCSELs 

We  have  recently  shown  that  discrimination  between  plasma  and  lattice  heating  is  important  in 
explaining  the  switch-off  of  VCSELs  under  increased  pumping  and  we  have  predicted  a  temperature 
dependent  hysteresis  effect  in  these  lasers  which  has  been  recently  observed  experimentally.  The 
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model  starts  from  the  microscopic  many-body  theory  and  builds  in,  in  an  incremental  manner,  the 
relevant  microscopic  effects  once  their  contribution  becomes  evident.  This  approach  lies  midway  in 
the  hierarchical  chain  alluded  to  above. 

•  Polarization  Patterns  in  Wide  Aperture  3-Level  and  Semiconductor  Lasers 

The  vector  character  of  the  laser  output  field  has  been  ignored  in  the  study  of  wide  aperture  lasers 
and  amplifiers.  It  is  weU  known  experimentally  that  wide  VCSEL  structures  show  polarization 
discrimination  in  the  transverse  plane.  Our  approach  is  to  develop  the  vector  order  parameter 
equation  description  for  basic  3-  and  4-  level  laser  models  with  the  latter  approximating  the  band 
structure  of  semiconductor  lasers.  The  additional  polarization  degree  of  freedom  considerably  extends 
the  available  range  of  experimentally  observable  behaviors.  The  combination  of  spatial  patterns  and 
polarization  domains  offers  exciting  opportunities  for  nonlinear  polarization  switching  and  control  of 
patterns  in  spatially  extended  lasers. 

•  ID  -  3D  Vector  Maxwell  Solvers 

Ultrashort  intense  light  pulses  propagating  through  wavelength  scale  inhomogeneous  structures  such 
as  periodic  index  gratings,  sharp  subwavelength  dimensional  interfaces,  Quantum  WeU  stacks.  Quan¬ 
tum  wire  structures,  semiconductor  micro-cavites  etc  aU  require  a  fuU  resolution  of  the  fuU  optical 
carrier  wave.  Possible  nonlinear  effects  associated  with  such  propagation  include  carrier  and  envelope 
shocking  of  intense  femtosecond  optical  pulses,  novel  nonlinear  localization  phenomena,  backseat ter- 
ing  from  nonUnearly  induced  index  inhomogeneities  etc.  Many  of  the  ID  solvers  can  be  optimized  to 
run  in  real  time  on  our  new  Power  ChaUenge  machine  offering  a  nonUnear  optical  CAD-like  environ¬ 
ment.  The  goal  of  this  project  is  to  develop  interactive  graphical  tools  which  maximize  efficiency  in 
transitioning  new  nonUnear  optical  effects  to  technology  appUcations. 

•  Dynamic  NonUnear  Optical  Skin  Effect 

This  project  utiUzed  the  vector  MaxweU  solvers  to  study  a  novel  nonUnear  reflection  phenomenon 
at  an  interface.  Our  prediction  is  that  an  intense  femtosecond  optical  pulse,  incident  on  a  nonUnear 
interface,  can  induce  a  moving  sub-wavelength  absorption  front  which  acts  as  a  curved  reflector. 
Direct  experimental  evidence  for  this  phenomenon  has  recently  been  provided  by  our  coUeagues  at 
the  Optical  Sciences  Center. 

•  Large  Scale  Computational  NonUnear  Optics 

Our  group  has  been  taking  the  lead  in  promoting  this  emerging  field  of  computational  science. 
The  recent  DURIP  award  of  an  SGI  Power  ChaUenge  L  plus  5  Indigo  graphics  workstations  has 
considerably  enhanced  our  intermediate  scale  supercomputing  and  real  time  graphical  visuaUzation 
capabiUties.  Our  two  main  large  scale  computational  projects  involve: 

(1)  developing  a  3D  +  time  vector  MaxweU  solver  on  the  CM5  at  AHPCRC,  Minnesota,  for  nonUnear 
optics  problems  of  relevance  to  the  present  project  and  to  a  companion  AFOSR  contract  No.  F49620- 
94-1-0051  “3D  CoUapse  Phenomena  in  Dispersive  NonUnear  Media:  A  Critique  of  Envelope  Models.” 

and 

(2)  developing  a  MaxweU-Semiconductor  Bloch  code  incorporating  fuU  many-body  interactions  on 
the  SP2  machine  on  Maui 

The  former  code  is  highly  paraUel  and  takes  advantage  of  these  features  on  the  CMS  whereas  the 
course-grained  paraUeUsm  of  the  latter  is  ideaUy  suited  for  the  SP2. 
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•  Relaxation  Processes  in  Semiconductor  Light  Amplifiers 

This  project  extends  the  description  of  the  nonequilibrium  ultrafast  relaxation  and  thermalization 
of  carriers  beyond  the  basic  quantum  Boltzmann  description.  The  idea  is  to  include  higher  order 
correlations  neglected  in  the  Kadanoff-Baym  theory.  As  semiconductor  lasers  are  examples  of  sys¬ 
tems  which  deviate  considerably  from  equilibrium,  the  role  of  higher  order  correlations  needs  to  be 
evaluated  in  this  context.  The  work  combines  the  method  of  nonlinear  optical  susceptibility  with  the 
multiple  scales  method  in  order  to  compute  the  next  order  correction  term  to  the  quantum  Boltzmann 
equation. 

•  Multimode  Laser  Dynamics 

Edge  emitting  semiconductor  lasers  typically  accommodate  20-30  longitudinal  modes  under  the  gain 
bandwidth.  Multi-longitudinal  mode  oscillations  are  also  observed  in  broad  area  semiconductor  lasers 
immediately  beyond  lasing  threshold.  This  project  explores  multi-longitudinal  mode  oscillations  in 
a  simpler  2-level  laser  with  and  without  external  mirror  feedback.  The  goal  is  to  use  our  experience 
gained  from  this  simpler  system  to  build  a  realistic  multi-longitudinal  mode  wide  aperture  semicon¬ 
ductor  laser  model.  The  rate  equation  approach  described  above  also  shows  multi-longitudinal  mode 
oscillations  in  a  broad  area  laser  but  suffers  from  the  fact  that  we  must  include  an  artificial  diffusion. 
The  present  approach  includes  the  fast  polarization  dynamics  and  the  problem  is  extremely  stiff 
mathematically, 

•  Circularly  Symmetric  Distributed  Feedback  Semiconductor  Lasers 

New  generations  of  circularly  symmetric  distributed  feedback  lasers  (CSDFB)  are  being  investigated 
at  a  preliminary  stage  using  a  standard  beam  propagation  algorithm.  These  lasers  offer  the  advan¬ 
tages  of  circularly  symmetric  output,  low- divergence  output  beam  and  the  ability  to  phase-lock  a  2D 
array. 

D*  RESEARCH  FACULTY  AND  GRADUATE  STUDENTS 

Principal  Investigator: 
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E.  TECHNICAL  REPORTS 


Femtosecond  Pulse  Propagation  in  Semiconductor  Amplifiers  and  Lasers 

We  have  performed  a  variety  of  computations  studying  the  behavior  of  femtosecond  pulses  in  semiconductor 
amplifiers.  We  employ  a  many-body  theory  describing  the  interaction  of  light  with  a  bulk  semiconductor 
material.  The  behavior  of  an  optical  pulse  with  full  width  at  half  maximum  (FWHM)  of  the  order  100 
-  500  femtoseconds,  propagating  over  a  long  distance  in  a  semiconductor  amplifier  has  been  studied  in 
both  strongly  confined  and  broad  area  planar  geometries.  The  role  of  the  various  physically  important 
interactions  on  the  pulse  evolution  along  the  amplifier  have  been  investigated. 

Our  initial  studies  reported  on  here  are  confined  to  a  strong  femtosecond  optical  pulse  interacting  with 
a  highly  excited  bulk  GaAs  semiconductor.  In  this  regime  the  equations  used  to  model  the  interaction 
describe  the  evolution  of  the  electric  field  envelope  E  in  the  traveling  frame  of  reference  ^  =  z  —  tcin 


BE 


koV 


(1) 


the  momentum  resolved  polarization  Pj  and  electron/hole  distribution  functions,  which  evolve  as 


dPg 

dt 


-i  {Ag  -  Wo)  Pg  -  iilginl  +  nj  -  1)  + 


dPg 

dt 


coll 


iP^Ctq  —  iPq^^  +  + 


coll 


(2) 

(3) 


with  the  renormalized  Rabi  frequency  Qg  =  T'j-g'P?'-  Here  and  A,  =  —  ^  J2g' 

denote  the  distribution  functions  for  electrons  or  holes  and  the  energy  dispersion  including  the  band  gap 
renormalization,  respectively.  The  Coulomb  potential,  V^,  is  treated  in  a  quasi-statical  screening  model. 
The  Rabi  frequency  Q  ==  •  E/ti  is  determined  by  the  dipole  matrix  element  dev  sind  the  amplitude  of 

the  external  electrical  field  E  = 

We  are  considering  a  broad  area  semiconductor  device.  The  lateral  “?/”  direction  has  strong  confinement 
due  to  the  layered  structure  of  the  amplifier.  This  has  the  effect  of  multiplying  the  right  hand  side  of 
equation  1  by  the  confinement  factor  F,  and  replacing  the  Rabi  frequency  Q  by  Ffi.  If  the  wave  is  also 
strongly  guided  in  the  transverse  direction,  a  correspondingly  smaller  confinement  factor  must  be  used, 
and  the  diffraction  term  V±E  may  be  ignored,  it  just  produces  an  effective  shift  in  Uq- 


Anomalous  absorption 

R.  Indik,  A.  Knorr,  R.  Binder,  J.V.  Moloney,  S.W,  Koch  and  W.W.  Chow 


Before  we  can  use  our  model  equations,  we  must  choose  some  specific  form  for  the  collision  terms.  To  start 
with  we  approximate  the  collision  terms  using  a  rate  equation  (no  memory)  with  uniform  time  constants 
=  (60  fs)-^ 


^Icoll  = 

^i„  = 

dt 


_2l  +  7lp 
2  « 


(4) 

(5) 


The  terms  T),  are  Fermi  distributions  with  chemical  potentials  and  temperature  T  chosen  so 

that  they  produce  the  same  carrier  densities  iV®/'*  =  dq  and  total  kinetic  energy  as  the  distributions 

n®/*.  This  form  for  the  collision  term  corresponds  to  considering  carrier-carrier  collisions,  all  of  which 
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preserve  total  kinetic  energy  and  the  numbers  of  electrons  and  of  holes.  The  Fermi  distributions  are 
assumed  to  have  a  common  temperature,  and  individual  chemical  potentials.  In  aU  of  our  computations, 
we  assume  the  polarizations  and  the  densities  are  isotropic  in  the  momenta  q. 

One  of  the  more  surprising  i)redictions  of  the  MSB  equations  is  that  if  one  increases  the  intensity  of  the  light 
pulse  sufficiently,  the  amplifier  will  become  an  absorber.  [1].  This  is  different  from  the  behavior  of  a  simple 
two  level  system,  where  one  can  saturate  the  gain,  but  can  never  convert  it  into  absorption.  The  source  of 
the  anomalous  absorption  is  the  coupling  of  the  light  field  with  the  electrons  in  higher  momentum  states. 
For  higher  momenta,  the  material  is  not  inverted,  i.e.  light  at  corresponding  (higher)  frequencies  would  be 
absorbed.  Even  light  at  the  lower  frequencies  is  absorbed  slightly.  Ordinarily,  that  absorption  is  more  than 
counterbalanced  by  the  gain  that  the  lower  frequency  light  experiences  as  a  result  of  its  interaction  with 
the  inverted  lower  momentum  states  that  are  resonantly  coupled.  However,  when  the  intensity  of  the  light 
is  sufficient  to  saturate  the  gain,  the  absorptive  coupling  to  the  higher  momentum  states  may  be  revealed. 
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Figure  1:  Evolution  of  intensity  (left)  and  density  (right)  during  propagation. 


H  one  assumes  for  simplicity  that  the  pulse  doesn’t  change  its  shape  as  it  propagates,  we  see  that  this  model 
predicts  that  there  is  an  intensity  Jq  such  that  whatever  the  initial  intensity,  after  sufficient  propagation  the 
intensity  wiU  return  to  Jq-  is  the  intensity  value  for  at  which  the  semiconductor  response  switches  from 
a  net  gain  I  <  Iq  to  a  net  absorption  if  I  >  Iq.  Thus  it  is  natural  to  ask:  what  is  the  long  term  behavior 
of  a  pulse  that  propagates  in  a  semiconductor  amplifier?  We  considered  this  question  [2]  and  found  yet 
another  surprising  prediction  of  this  model.  Figure  1  shows  the  result  of  propagating  a  relatively  weak 
pulse  in  such  an  amplifier  over  a  long  distance.  Initially  the  pulse  is  amplified  linearly,  and  we  can  see  slight 
effects  of  gain  dispersion  and  group  velocity  dispersion  in  the  reshaping  of  the  pulse  and  its  spectrum.  This 
continues  until  the  amplifier  is  saturated.  This  saturation  is  a  consequence  of  plasma  heating,  spectral 
hole  burning  and  total  carrier  density  depletion.  Carrier-carrier  collisions  ‘‘refiU”  the  spectral  hole,  but 
the  model  in  equations  11  and  12  have  no  way  for  the  plasma  to  cool,  except  through  stimulated  emission 
and  absorption.  Once  the  plasma  is  heated,  the  carriers  are  distributed  across  more  momentum  states, 
and  fewer  of  those  states  are  inverted.  Plasma  heating  is  the  dominant  source  of  gain  saturation.  As  is 
nearly  always  the  case  with  a  saturable  amplifier,  the  leading  portion  of  the  pulse  steepens,  and  the  peak 
of  the  pulse  actually  moves  faster  than  the  speed  of  light  in  the  material.  In  addition,  the  refractive  index 
dispersion  is  such  that  the  redder  frequencies  travel  faster  than  the  blue.  Thus  these  frequencies  end  up 
closer  to  the  leading  edge  of  the  pulse,  and  experience  more  gain.  Thus  we  see  the  center  of  the  spectrum 
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of  the  pulse  shifting  below  the  peak  of  the  gain  spectrum.  Self-phase  modulation  also  contributes  to  this 
eflFect.  The  reshaping  of  the  pulse,  its  splitting  into  a  very  short  very  intense  leading  portion  and  much 
weaker  and  broader  tail  is  harder  to  understand.  We  can  see  by  comparing  the  shape  of  the  intensity  of  the 
light  pulse  with  the  total  carrier  density,  that  we  have  entered  a  regime  of  adiabatic  following.  Observation 
of  the  q  resolved  density  and  polarization  shows  that  much  of  that  adiabatic  following  is  due  to  coupling  to 
very  high  momentum  states.  The  robustness  of  this  phenomenon  was  checked  by  varying  the  amplitude, 
shape  and  frequency  of  the  initial  pulse.  Without  exception,  if  we  propagated  far  enough,  the  pulse  entered 
a  regime  of  adiabatic  following.  This  observation  was  intriguing  as  we  felt  that  the  presence  of  a  very  short 
intense  pulse  that  forms  spontaneously  and  robustly  could  be  useful  in  applications. 


Figure  2:  Intensity  (left)  and  Density  (right)  as  functions  of  t  and  x. 

The  transverse  confinement  was  then  removed  to  study  the  propagation  of  pulse  with  a  Gaussian  transverse 
cross-section  over  a  distance  of  1mm  (length  leading  to  a  linear  gain  for  the  intensity  of  e^°).  It  is  now 
observed  that  the  pulse  produces  its  own  transverse  waveguide.  The  intense  central  portion  depletes  the 
carrier  density.  This  lower  carrier  density  induces  a  higher  refractive  index,  and  therefore  a  focusing  effect. 
Figure  2  shows  the  amplitude  and  the  carrier  density  as  a  function  of  time  and  the  transverse  dimension 
after  the  beam  has  propagated  1mm.  Figure  3  shows  a  comparison  of  the  amplitude  in  the  center  of  the 
beam  with  transverse  structure  and  the  amplitude  of  a  beam  with  strong  transverse  guiding. 


lime  (ps) 

Figure  3:  Comparison  of  intensity  of  planewave  (solid)  and  center  value  of  intensity  for  a 
Gaussian-beam  after  propagating  1  mm. 
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We  have  also  explored  the  behavior  of  the  amplifier  when  two  strong  pulses  at  different  frequencies  are 
compete  for  gain  [3].  In  tha.t  case  we  see  the  effects  of  plasma  heating,  four  wave  mixing  and  spectral  hole 
burning.  Figure  4  shows  the  evolution  of  the  spectra  of  such  interacting  pulses. 
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Figure  4:  Evolution  of  the  spectrum  for  pairs  of  coincident  pulses  with  different  frequencies. 
Dashed  lines  show  evolution  of  individual  pulses. 


Plasma  cooling 

R.  Indik,  R.  Binder,  M.  Mlejnek  and  J.V,  Moloney 

One  of  the  more  important  features  governing  the  evolution  of  the  propagating  pulse  above  is  that  the 
plasma  cannot  cool  other  than  through  interaction  with  the  light  field  itself.  On  the  whole,  the  interaction 
with  the  light  field  would  tend  to  heat  the  plasma,  since  the  absorption  tends  to  be  into  high  momentum 
states,  (creating  fast  carriers),  and  the  emission  removes  low  momentum  carriers.  The  collision  terms, 
however,  should  include  terms  we  initially  neglected  that  tend  to  cool  the  plasma.  The  carriers  collide 
not  only  with  each  other,  but  also  with  optical  phonons.  The  optical  phonons  can  take  energy  from  the 
carriers  and  couple  it  to  the  lattice.  This  process  is  slower  than  carrier- carrier  scattering,  but  is  stiU  fairly 
rapid.  By  using  a  quasi-equilibrium  assumption  one  can  estimate  the  cooling  rate  for  the  carriers.  Once 
again,  we  add  this  term  to  our  equation  using  a  rate  equation  approximation. 


^IcoU 


dn^’' 

dt 


7r +7, 


ho 


-if  [nf  -  -  if^o  [nf  - 


(6) 

(7) 


Here  we  have  introduced  two  new  rates  and  two  new  Fermi  distributions.  The  new  Fermi  distributions 
are  chosen  to  preserve  the  total  carrier  density,  but  have  a  fixed  lattice  temperature  Tj^,  On 
much  longer  time  scales,  one  must  consider  that  the  lattice  itself  will  be  heated,  and  thus  dynamics  for  Tjr 
could  become  important  [4].  This  much  slower  thermal  process  is  not  relevant  to  our  study.  As  expected, 
the  effect  of  including  the  cooling  of  the  plasma  is  to  make  the  saturation  of  the  gain  more  difficult.  As  the 
plasma  cools,  carriers  move  from  high  momentum  states  to  lower  ones  and  once  again  become  available 
for  stimulated  emission.  In  thin  slabs,  when  the  pulses  are  short  enough,  there  is  little  effect.  Table  1 
shows  how  the  intensity  (presented  in  terms  of  the  corresponding  Rabi  frequency  Do)  at  which  the  gain  is 
saturated  varies  with  pulse  length  with  and  without  cooling. 
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Pulse  length  (FWHM) 

Qo  (iio  cooling) 

fio  (with  cooling) 

100  fs 

21.89  mev 

22.85  mev 

200  fs 

16.39  mev 

18.81  mev 

400  fs 

12.03  mev 

16.49  mev 

800  fs 

8.67  mev 

15.78  mev 

1600  fs 

6.18  mev 

17.78  mev 

Table  1:  J2o  with  and  without  cooling  as  a  function  of  the  pulse  length. 

For  longer  propagation  lengths,  we  see  that  cooling  has  a  dramatic  eflPect  even  for  pulses  that  are  initially 
very  short.  Where  previously  the  gain  was  depleted  by  the  leading  portion  of  the  pulse,  we  now  see  that 
the  cooled  carriers  produce  additional  gain  on  the  trailing  portion  of  the  pulse.  The  leading  edge  of  the 
pulse  is  not  greatly  changed  (at  least  qualitatively),  but  now  the  vast  bulk  of  the  energy  in  the  pulse  is 
on  the  trailing  portion.  To  take  advantage  of  the  very  sharp  leading  portion  of  the  pulse,  it  would  have 
to  be  separated  from  its  much  larger  tail.  Figure  5  shows  the  time  resolved  intensity  of  the  pulse  after 
propagation  roughly  2mm.  Figure  6  shows  the  transverse  structure  formed  by  a  150  femtosecond  pulse  with 
a  transverse  Gaussian  profile  after  propagating  1mm.  Note  that  the  self  guiding  is  stiU  quite  pronounced. 
We  have  also  looked  at  the  effect  of  varying  the  shape  of  the  initial  pulse  on  the  output  of  the  long  amplifier. 
In  contrast  to  the  case  without  optical  phonon  collisions,  we  find  that  the  shape  of  the  output  varies  with 
the  shape  of  the  input.  The  “wiggles”  that  foUow  the  sharp  initial  peak  in  the  output  are  more  pronounced 
for  initial  shapes  that  fall  off  more  quickly. 


time  (ps) 

Figure  5:  Intensity  of  pulse  after  propagating  2mm.  Cooling  effects  have  been  included. 
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Memory  effects 

R.  Indik,  R.  Binder,  Moloney,  A.  Knorr,  S.W.  Koch 

The  rate  equation  form  that  we  choose  for  the  collision  term  in  the  decay  of  the  polarization  is  a  result 
of  making  a  simplifying  assumption  about  carrier-carrier  scattering.  Namely,  we  assume  that  the  colli¬ 
sion  process  has  no  memory.  There  is  significant  experimental  evidence  to  suggest  the  contrary.  In  the 
low  density  case,  it  is  knov;n  that  the  absorption  falls  of  exponentially  with  detuning  below  the  exciton 
frequency.  This  “Urbach”  tail  has  been  modeled  using  a  Gaussian  or  a  sech  lineshape  rather  than  the 
Lorentzian  profile  consistent  with  a  fixed  polarization  damping  rate.  Even  at  higher  densities,  we  suspect 
that  the  Lorentzian  line  shape  (corresponding  to  no  memory)  must  be  incorrect.  Below  the  (renormalized) 
bandgap,  there  is  very  little  absorption.  However,  in  a  bulk  semiconductor,  use  of  a  Lorentzian  lineshape 
leads  to  a  gain  spectrum  with  an  absorptive  tail  that  extends  well  below  the  bandgap.  If  a  sech  lineshape 
is  used,  this  feature  disappears.  If  we  examine  the  argument  that  indicates  how  the  high  momentum  states 
together  with  the  quadratic  density  of  states  cumulatively  produce  significant  absorption  for  frequencies 
that  ordinarily  would  see  net  gain,  we  can  see  that  we  need  only  assume  that  the  lineshape  falls  oflF  faster 
than  1/^^  to  remove  the  tails.  Thus  a  decrease  like  is  enough  to  greatly  diminish  the  influence  of  the 
high  momentum  states.  There  is  some  literature,  proceeding  from  first  principles  using  quantum  Boltz¬ 
mann  equations,  that  calculates  lineshapes  with  asymptotics  as  described  above.  Such  models  are  quite 
significantly  more  complex  than  those  we  have  been  using,  and  we  have  captured  the  essential  features  (i.e. 
the  aymptotics)  using  a  phenomenological  approach. 

It  turns  out  that  if  we  take  for  the  collision  term  of  the  polarization 


dP, 

dt 


coll 


/  Qgi*  -  t')Pq(t') 

Jo 


(8) 


where 

the  lineshape  will  have  the  desired  quartic  faU  off  with  detuning,  and  that  if  we  restrict  to  such  simple 
diagonal  memory  terms  this  is  essentially  unique.  We  caU  this  a  two  pole  model,  as  in  Fourier  space  the 
response  function  corresponding  to  the  diagonal  linear  terms  in  the  Pq  can  be  written  as  a  rational  function 
with  two  poles  in  the  complex  upper  half  plane.  One  very  convenient  feature  of  such  a  form,  is  that  it  is 
easily  converted  back  into  a  system  of  coupled  ODEs.  We  replace  the  P  equation  with  a  pair  of  equations, 


Pg  =  WgPg-'flgQg  +  iOg{nl  +  n^^-l)  (9) 

Qi  =  {-llq-l2q+i<^q)Qq-i-l2qPq  (10) 


Assuming  a  two  pole  form  for  the  response,  is  equivalent  to  replacing  the  Markovian  assumption  that  leads 
to  the  Lorentzian  line  shape  with  a  two-stage  Markovian  assumption. 


detuning  _(mev)  detuning 

Figure  7 :  Comparison  of  gain  spectra  using  Lorentzian  (dash-dot ),  sech  (solid),  and  two  pole 
(dashed)  models  (left),  and  lineshapes  (right) 
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Figure  7  (right)  shows  a  comparison  of  the  gain  computed  using  the  Lorentzian,  several  choices  of  parameter 
values  of  the  2-pole  model,  and  a  sech  lineshape.  Our  preferred  set  of  parameters  is  chosen  to  make  the 
2-pole  lineshape  approximate  the  sech  lineshape.  The  Lorentzian  tail  is  removed  by  both  the  sech  and 
the  2  pole  models.  In  fact  there  is  little  to  distinguish  them  when  looking  at  the  gain,  though  both  differ 
significantly  from  the  Lorentzian.  Figure  8  (left)  compares  the  response  of  the  total  density  to  a  very  strong 
pulse.  Cooling  effects  have  not  been  included  as  they  are  not  important  for  such  short  pulses  in  a  thin  slab. 
We  see  that  the  anomalous  absorption  is  absent.  This  should  be  understood  in  terms  of  the  decreased 
absorptive  coupling  to  the  high  momentum  states.  Instead  we  see  a  new  feature.  The  density  seems  to 
follow  the  field  intensity  up  and  down,  temporarily  increasing  over  its  initial  value,  while  the  strong  pulse  is 
present.  This  behavior  can  be  understood  from  the  behavior  of  the  real  part  of  the  susceptibility  x'(‘^)-  By 
construction,  the  imaginary  part  of  the  susceptibility  decays  rapidly  with  detuning,  but  applying  Kramers- 
Kronig  shows  us  that  the  real  part  stiU  decays  like  the  inverse  of  the  detuning.  This  is  the  coupling  that 
leads  to  adiabatic  following.  This  is  just  as  strong  as  before,  but  now  the  absorptive  coupling  is  much 
weaker,  and  the  adiabatic  following  is  dominant.  If  we  increase  the  intensity  of  the  pulse  a  great  deal,  we 
can  once  again  see  anomalous  absorption,  however,  the  intensities  are  now  quite  huge,  and  the  adiabatic 
following  is  now  very  pronounced  (Figure  8  right).  We  have  also  calculated  the  result  of  propagation 
without  cooling,  but  using  the  2-pole  model.  It  is  interesting  that  these  memory  effects,  even  without 
cooling,  provide  gain  behind  the  leading  edge  of  the  pulse.  Figure  9  shows  the  corresponding  results  for 
propagation  in  a  broad  area  amplifier  using  the  2-pole  model  with  cooling  included.  The  similarity  to  the 
plane  wave  propagation  due  to  the  presence  of  the  self  trapping  is  stiU  quite  marked. 
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Figure  8:  (left)  Comparison  of  response  of  the  total  density  to  a  very  strong  using  Lorentzian  (upper  solid), 
sech  (lower  solid),  and  and  various  two  pole  models  (dashed),  (right)  response  of  the  two  pole  model  to  a 
pulse  9  times  as  intense. 


Conclusions 

Our  preliminary  studies  of  femtosecond  pulse  propagation  in  a  bulk  semiconductor  amplifier  have  revealed 
a  subtle  and  complex  interplay  between  various  physical  mechanisms  which  can  be  ascribed  to  many- 
body  effects.  The  significant  contribution  of  these  phenomena  for  femtosecond  probes  even  stretches  the 
bounds  of  the  many-body  theory  as  currently  developed  and  it  is  highly  unlikely  that  simple  parameterized 
rate  equation  models  can  provide  any  significant  insight.  In  particular,  we  have  studied  systematically 
the  role  of  carrier-carrier  scattering,  carrier-phonon  scattering  and  memory  effects  on  femtosecond  pulse 
propagation  in  narrow  stripe  and  broad  area  lasers.  Plasma  cooling  and  memory  effects  are  shown  to  lead 
to  a  pronounced  amplification  of  the  trailing  edge  of  the  pulse  and  hence  severe  pulse  distortion  both  in 
single  stripe  and  broad  area  structures.  In  addition  these  effects  accentuate  the  self-induced  transient  index 
guide  formed  in  a  broad  ar<5a  structure. 

It  is  expected  that  these  many-body  ingredients  will  allow  us  to  evaluate  the  ultimate  modulation  response 
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of  mode-locked  semiconductor  lasers  and  unearth  the  source  of  the  observed  strong  frequency  chirp  in 
these  devices.  Computationally,  this  class  of  problem  poses  significant  challenges  especially  in  modeling 
broad  area  amplifier  structures  or  lasers.  On  the  other  hand,  it  will  allow  us  to  build  more  manageable 
and  realistic  semi-phenomenological  amplifier  and  laser  models  based  on  the  insight  gained  from  the  more 
comprehensive  many-body  approach. 


Figure  9:  Transverse  structure  of  intensity  (left)  and  density  (right)  after  1mm  when  cooling 
and  memory  effects  are  included. 
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Microscopically  Computed  Quantum  Well  Gain  Spectra 

S.  Hughes,  A.  Girndt,  C.Z.  Ning,  S.W.  Koch,  R.  Indik,  R.  Binder  and  J.V.  Moloney 

An  open  question  that  remains  in  the  study  of  many-body  interactions  in  semiconductor  media  is  how  to 
effect  a  closure  of  the  infinite  hierarchy  of  equations  describing  the  nonlinear  optical  response.  A  significant 
discrepancy  existed  between  experimentally  measured  and  microscopically  computed  quantum  well  gain 
spectra  until  a  few  months  ago.  We  have  recently  identified  this  discrepancy  with  an  inconsistent  treatment 
of  carrier  collision  effects,  at  the  level  of  quantum  kinetic  theory  in  the  Markovian  limit.  Earlier  comparisons 
between  bulk  spectra  where  the  coUison  terms  in  the  Semiconductor  Bloch  equations  (SBE)  were  treated 
in  a  rate  equation  approximation,  and  experiment  proved  to  be  excellent.  An  important  conclusion  to  be 
drawn  from  this  study  is  that  femtosecond  carrier-carrier  and  carrier-phonon  scattering  events,  profoundly 
affect  the  shape  of  the  CW  semiconductor  optical  response  function. 

In  the  microscopic  treatment  of  the  semiconductor  optical  response,  the  momentum  resolved  polarization 
Pk  and  electron/hole  distribution  functions,  evolve  as 


dt 

dn]!'^ 

dt 


-i  (Ak  -  Wo)  Pk  -  iilkinl  +  «*-!)  + 
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with  the  renormalized  Rabi  frequency  fit  =  f  +  ^  Ef  Vk-k'Pk'-  Here  nli'‘  and  A*  =  e*  -  i  Et-  14-J:'(^^^  + 
njj,)  denote  the  distribution  functions  for  electrons  or  holes  and  the  energy  dispersion  including  the  band 
gap  renormalization,  respectively.  The  Coulomb  potential,  14,  is  treated  in  a  quasi-statical  screening  model. 
The  Rabi  frequency  fi  =  d^v  •  E/h  is  determined  by  the  dipole  matrix  element  d^v  and  the  amplitude  of 
the  external  electrical  field  E  = 

The  issues  of  anomalous  absorption,  memory  effects  etc,  discussed  above  can  be  traced  to  how  one  computes 
the  collision  terms  in  the  polarization  equation.  When  calculating  semiconductor  gain  spectra,  one  often 
approximates  collision  effects  with  a  constant  effective  dephasing  rate,  which  leads  to  anomalous  absorption 
below  the  renormalized  bandgap.  Moreover,  experimentally  measured  quantum- weU  gain  spectra  look  bulk- 
Uke  whereas  computed  gain  spectra  show  a  sharp  rise  near  the  band  edge[l].  A  more  extended  treatment 
of  collision  effects  involves  extending  the  screened  Hartree-Fock  treatment  to  include  contributions  from 
the  next  higher  order  correlations  involving  the  Coulomb  potential[2,  3].  The  resulting  equation  for  the 
microscopic  polarization  due  to  an  electron- hole  pair,  pk,  has  the  form 

Q 

—Pk  =  -iiVkPk  -  i^k{n^,k  +  nh,k  -  1)  -  TkPk  +  '^Tk,,Pk+g.  (13) 

The  first  two  terms  on  the  right  hand  side  describe  the  oscillation  of  the  polarization  at  the  transition 
frequency,  a?*,  and  the  stimulated  emission  and  absorption  processes.  The  screened  Hartree-Fock  contri¬ 
butions  lead  to  a  bandgap  renormalization,  resulting  in  the  above  transition  energy  and  the  renormalized 
Rabi  frequency  above.  Carrier-carrier  collisions  give  rise  to  the  last  two  terms.  The  third  term  is  a  diagonal 
contribution,  with 


27r 


r*=  E  g^(^a,k  +  €b,k-€a,k+q  —  €h,k~g)'X[na,k+qi^-nt^k)nb,k-q  +  i'i-  —  na,k+g)ni,k{'i-  —  ni,  *_«)] 

- -  I.  I _ L  y, 


a:=e,h  b=e,h  q 


(14) 

The  fourth  term  is  a  nondiagonal  contribution,  E,  which  couples  the  polarizations  with  different 

A:’s.  For  this  contribution. 
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^  ^  y  ^  ^a,fc+3  — })^[(1  ^o,t)(l  ^ft,<;)^J,ife— }"t"^o,fc^d,fc(l  ^b,k—q}^  (^^) 

a=e,h  b=e,h  k  ^ 

The  dramatic  modification  of  the  quantum  well  gain  spectrum  when  the  constant  dephasing  assumption 
is  removed,  is  illustrated  in  Figure  1(a)  where  the  gain  spectra  with  (dashed)  and  without  (solid)  this 
assumption  are  compared  for  three  values  of  the  total  carrier  density.  Note  that  the  sharp  rise  on  the 
negative  detuning  side  is  removed  and  the  peak  gain  values  are  significantly  larger.  Figure  1(b)  is  a 
comparison  of  the  new  gain  spectra  with  experimental  measurements  of  D.  Bossert  at  PL/LIDA.  The  solid 
squares  are  experimental  measurements  and  the  agreement  is  seen  to  be  excellent.  This  latter  picture  is 
reproduced  in  Color  Plate  I  below. 

The  computation  of  the  full  gain  spectra  requires  supercomputing  resources  but  the  payoff  is  huge.  Now, 
for  the  first  time,  we  have  confidence  that  we  can  design  a  laser  structure  from  the  bottom  up  and  predict  its 
performance  prior  to  fabrication.  The  gain  spectra  have  recently  been  computed  for  other  QW  structures 
and  verified  experimentally  by  a  different  group[4]. 
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Universal  order  parameter  equation  description  of  laser  patterns 
J.  Lega,  Q.  Feng,  J.V.  Moloney  and  A.C.  Newell 

This  project  derives  the  complex  order  parameter  equations  for  a  laser  starting  with  the  physical  laser  equa¬ 
tions.  A  goal  is  to  evaluate  the  robustness  of  the  order  parameter  equation  model  relative  to  the  original 
physical  model  as  a  precursor  to  deriving  an  order  parameter  equation  description  for  a  broad  area  semi¬ 
conductor  laser.  The  latter  would  start  from  the  full  many-body  Maxwell- Semiconductor  Bloch  equations 
and  some  preliminary  progress  has  been  made  on  this  as  noted  below.  Our  prototypical  models  have  been 
the  2-level  homogeneously,  inhomogenously  broadened  and  Raman  (3-level)  lasers.  These  systems  sepa¬ 
rately  mirror  a  particular  physical  manifestation  of  the  known  behavior  of  the  semiconductor  laser  at  a  less 
complicated  level.  The  latter  is  assumed  to  have  a  2-band  structure  and  is  inhomogeneously  broadened  in 
electron  momentum  space.  In  a  semiconductor,  the  carrier  density  dependent  detuning  dynamically  alters 
the  bandgap  (bandgap  renormalization),  while  a  Raman  laser  exhibits  an  intensity  dependent  detuning 
(A.C.  Stark  effect).  By  studying  the  role  of  the  constituent  parts  in  simpler  mathematical  settings  we 
hope  to  glean  the  influence  of  analogous  physics  in  the  much  more  complicated  many-body  semiconductor 
problem. 

We  have  derived  the  complex  Ginzburg-Landau  (CGL)(  negative  detuning),  and  coupled  complex  NeweU- 
Whitehead-Segel  (CNWS)  (positive  detuning)  amplitude  equations  for  both  the  2-level  and  Raman  lasers 
[1].  Very  recently,  we  derived  a  complex  Swift-Hohenberg  (CSHE)  description[3]  which,  although  strictly 
valid  near  the  peak  gain  of  the  laser,  proves  to  be  very  robust.  The  CSH  model  agrees  with  the  numericaffy 
computed  instability  curves  for  the  2-level  laser  well  beyond  lasing  threshold.  The  CSH  has  been  extended 
to  derive  a  ’’generalized  rate  equation”  for  a  mathematically  stiff  laser  which  would  represent  a  lowest 
order  approximation  to  a  realistic  semiconductor  laser  model.  This  extended  system  is  given  by: 

-I-  - ~r(^  “  aV^)^V’  “ 

<7  -h  1 

Tit  =  -bn  -1- 1^1^ 

This  coupled  system  has  an  obvious  physical  interpretation.  The  linear  growth  rate  (gain)  is  represented  by 
/i,  the  Hopf  bifurcation  frequency  (frequency  pulling)  by  ft  ,  the  diffraction  parameter  is  a  ,  cr  is  the  cavity 
loss,  and  b  is  the  normalized  inversion  decay  rate.  The  diffusion  term  (third  term  on  the  right)  acts  as 
a  gain  discrimination,  favoring  the  emission  of  an  off-axis  far-fleld  output  at  critical  transfer  wavenumber 
kc  =  by  damping  modes  at  either  side  of  this  wavennumber.  Without  this  term,  spurious  high 

wavenumber  modes  sees  growth  and  look  like  numerical  grid  osciUationsjl]. 

The  mean  flow  term  given  by  the  second  equation  is  important  in  contributing  to  the  robustness  of  the 
model.  With  this  equation  appended,  we  find  that  the  various  phase  and  amplitude  instability  boundaries 
of  this  ’’rate  equation”  agree  quantitatively  with  the  full  MaxweU-Bloch  2-level  laser  equations  up  to  two 
to  three  times  threshold[4].  It  is  easy  to  demonstrate  in  fact  that  the  single  CSHE  is  insensitive  to  the 
stiffness  of  the  system  (parameter  6)  and  that  the  additional  equation  is  therefore  essential.  We  have  also 
shown  that  the  usual  amplitude  equations  derived  in  [2]  for  negative  (complex  Ginzburg-Landau  (CGL)) 
and  positive  (complex  Newell- Whitehead- Segel  (CNWS))  detuning  from  peak  gain  are  in  precise  agreement 
with  those  derived  directly  from  the  above  coupled  CSHE  model. 

An  important  conclusion  from  this  study  is  that  the  stiffness  inherent  in  semiconductor  lasers  causes  the 
Busse  balloon  (stable  lasing  emission  domain)  to  disappear,  implying  that  the  laser  is  intrinsically  unstable 
even  at  lasing  threshold.  This  is  consistent  with  experimental  observations  showing  strong  uncontrolled 
filamentation  in  semiconductor  lasers.  Mathematically,  this  instability  can  be  ascribed  to  a  strong  pertur¬ 
bation  by  the  eigenvalue  of  the  linearization,  associated  with  the  mean  flow  (n),  of  the  complex  eigenvalues 
associated  with  the  order  parameter  tp.  Although  this  eigenvalue  can  be  damped  at  A;  =  0,  its  real  part  can 
become  positive  at  finite  k. 
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We  have  also  investigated  the  influence  of  highly  absorbing  boundaries  on  pattern  selection  in  these  large 
aspect  ratio  lasers.  One  finds  that  the  boundary  forces  a  selected  wavenumber  on  the  pattern  in  the  bulk[5]. 
A  quantitative  understanding  of  these  issues  is  vital  to  understanding  how  to  control  or  usefuUy  manipulate 
filamentation  instabilities  in  broad  area  and  flared  high-power  semiconductor  lasers.  In  addition  we  know 
that  these  lasers  run  on  multi-longitudinal  modes  even  just  above  threshold.  The  state  of  the  art  of 
semiconductor  laser  modeling  in  the  general  research  community  is  still  based  on  the  Beam  Propagation 
Method  which  ignores  time  dependence  entirely.  Simply  extending  this  method  by  adding  time  dependence 
wiU  lead  to  the  spurious  oscillations  mentioned  above. 

This  work  is  provides  a  vital  input  to  our  modeling  of  realistic  wide  aperture  semiconductor  lasers  as 
described  below.  Proof  of  principle  methodology  can  be  developed  within  the  present  framework  utilizing  a 
mix  of  analytical  and  computational  approaches.  Extending  the  current  model  to  include  multi-longitudinal 
modes  is  a  high  priority. 
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Controlling  Optical  Turbulence 
D*  Hochheiser,  J.V*  Moloney  and  J.  Lega 

This  project  presents  a  robust  global  control  strategy  for  a  class  of  partial  diflferential  equations  exhibiting 
spatio-temporal  (iisor(ier[l].  The  control,  implemented  as  a  spatial  filter  with  delayed  feedback,  is  shown  to 
stabilize  and  steer  the  weaWy  turbulent  output  of  a  spatially-extended  system.  The  latter  is  described  by 
a  generalized  complex  Swift-Hohenberg  equation  [2]  which  is  used  as  a  generic  model  for  pattern  formation 
in  the  transverse  section  of  semiconductor  lasers. 

The  technique  we  propose  takes  advantage  of  the  fact  that  the  spatio-temporal  spectrum  of  the  turbulent 
output  often  satisfies,  to  a  good  approximation,  the  dispersion  relation  of  the  system.  Control  is  then 
achieved  by  filtering  this  spectral  output  about  a  desired  plane  wave  and  feeding  back  the  corresponding 
delayed  signal  into  the  system,  together  with  a  contribution  proportional  to  the  retarded  field.  This 
method  is  likely  to  be  applicable  to  any  partial  differential  equation  which  sustains  traveling  waves  and  is 
particularly  suited  for  optical  systems.  In  the  latter,  the  farfield  output  is  indeed  a  natural  spatial  Fourier 
transform  of  the  complex  envelope  of  the  electromagnetic  field  and  applying  a  temporal  Fourier  transform 
to  successive  readings  of  the  farfield  yields  the  desired  spatio-temporal  spectrum. 

The  control  scheme  has  been  illustrated  on  the  example  of  a  generalized  laser  Swift-Hohenberg  equation 
[2],  which  phenomenologically  describes  the  dynamics  of  wide  aperture  semiconductor  lasers.  The  Swift- 
Hohenberg  equation  is  a  well-known  generic  model  of  pattern  formation  in  extended  systems,  and  it  was 
shown  in  [2]  that  when  coupled  to  an  equation  for  the  population  inversion,  a  complex  SH  equation  gives  a 
good  description  of  2-level,  class  B  lasers.  Particular  effects  can  then  be  included  in  this  model  to  reproduce 
filament  at  ion,  a  typical  feature  of  wide  aperture  semiconductor  lasers.  The  latter  are  interesting  physical 
manifestations  of  spatially  extended  systems  showing  persistent  weakly  turbulent  behavior.  In  contrast 
to  wide  aperture  two-level  lasers,  they  display  strong  dynamic  filamentation  instabilities  immediately  at 
threshold.  Moreover,  their  large  gain  leads  to  strongly  amplified  spontaneous  emission  along  the  laser  axis 
and  hence,  to  a  persistent  noisy  background  behavior. 

Ideally,  a  robust  control  strategy  for  an  infinite  dimensional  nonlinear  dynamical  system  should  require  no 
a  priori  knowledge  of  the  solution  to  be  stabilized  and  should  take  into  account  the  fact  that  the  system 
may  have  no  stable  attracting  state.  The  method  we  propose  here  not  only  stabilizes  unstable  traveling 
waves  in  the  turbulent  regime  but  allows  one  to  select  and  angle  tune  (steer)  the  system  output  starting 
from  initial  noise  or  a  turbulent  state. 

The  most  important  distinction  between  the  semiconductor  and  other  laser  systems  is  the  marked  asymme¬ 
try  of  its  gain  and  refractive  index  spectra.  This  causes  a  very  strong  nonlinear  amplitude-phase  coupling 
in  the  field  which  leads  to  uncontrolled  dynamical  filamentation  in  the  laser  intensity  at  and  beyond 
threshold.  This  effect  is  modeled  by  a  nonlinear  coupling  between  the  electric  field  and  the  carrier  density 
through  a  coefficient  known  as  the  a-  factor.  By  introducing  a  similar  term  in  the  2-level  laser  complex 
Swift-Hohenberg  (CSH)  equation,  we  obtain  the  following  system: 

=  (T{r{x)  -  1)^  +  -  iaQ'ip 

+  <T{l  +  ia)ni;  (1) 

=  -fen+IV’P, 

where  a  is  negative  for  semiconductor  lasers. 

The  complex  order  parameter  'll;  is  the  scaled  envelope  of  the  electric  field  and  n  is  a  scaled  relative 
carrier  density.  The  latter  acts  as  a  mean-flow  and  has  a  profound  influence  in  destabilizing  the  system, 
leading  to  a  very  complicated  linear  growth  behavior  of  the  traveling  wave  solutions.  Here  a  is  the  scaled 
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cavity  loss  coefficient,  a  is  proportional  to  the  inverse  of  the  Fresnel  number  of  the  laser  and  measures 
the  characteristic  length  scale  in  the  transverse  dimension  relative  to  the  wavelength  of  light,  ft  is  the 
dimensionless  detuning  of  the  laser  frequency  from  the  gain  peak,  is  generally  the  two  dimensional 
Laplacian  although  we  restrict  our  study  here  to  one  transverse  dimension  x,  and  b  is  the  dimensionless 
ratio  of  the  carrier  recombination  to  polarization  dephasing  times  in  the  Semiconductor  Bloch  equations. 
The  external  pump  parameter  r(x)  is  the  scaled  external  current  applied  to  the  laser  and  the  x-dependence 
is  explicitly  displayed  in  order  to  emphasize  that  the  pumping  is  only  applied  over  a  finite  transverse  section 
of  the  laser.  Ontside  the  pumped  region,  the  passive  semiconductor  acts  as  a  very  strong  absorber.  In  the 
discussion  below  we  first  assume  that  the  pump  is  infinitely  extended  in  x,  in  order  to  take  advantage  of  the 
known  properties  of  solutions  to  the  CSH  equation.  We  then  show  explicitly  that  many  of  the  properties 
of  this  idealized  system  carry  over  to  the  realistic  finitely  pumped  cross-section. 

Equation  (1)  admits  traveling  wave  solutions  of  the  form  tj)  =  jRe*(*  n  =  R^/b  where  R?  =  b{r  —  1  — 
frequency  of  the  traveling  wave  of  wavenumber  k  is  given  by  (1  -|-  (7)w  =  — -f  ftcr  -|- 
aa[r  —  1  —  )^]]-  This  latter  expression  shows  that  the  nonlinear  amplitude- phase  coupling  makes  the 

frequency  dependent  on  the  distance  above  threshold  for  lasing. 

Optically,  as  the  nonlinear  dispersion  curve  {or  k  —  u  spectrum)  is  an  experimentally  accessible  quantity,  it 
is  natural  to  introduce  an  optical  feedback  scheme  which  filters,  in  k  and  u,  the  desired  traveling  wave  lying 
on  the  actual  dispersion  curve  (attractor).  The  idea  is  then  to  introduce  a  delayed  feedback  with  a  spatial 
filter  (consisting  of  a  lens  and  aperture  at  the  focal  point  of  the  lens).  Ideally  the  delay  (proportional 
to  w“^)  should  be  chosen  to  locate  the  desired  coordinate  on  the  k  —  u  plot.  If  this  scheme  works,  we 
simultaneously  have  achieved  stabilization  and  a  beam  steering  capability  by  simply  tuning  the  filter  along 
the  experimental  dispersion  curve.  Our  control  technique  is  then  the  following:  we  add  a  feedback  term  of 
the  form  — 7(V’(x,/)  —  '^(x,t  —  r))  to  the  above  equation  where 

'ip{x,t-T)  =  {F{k  -  ko)  T {ipix,  t  -  r))) , 

T  is  the  Fourier  transform  operator  and  jF  is  a  suitable  aperture.  This  represents  a  time-delayed,  spatially 
filtered  feedback  of  the  original  complex  field  at  the  output  facet. 

With  this  technique,  the  unstable  region,  which  for  a  <  0  corresponds  to  the  whole  domain  of  existence 
of  traveling  waves,  can  be  stabilized  over  the  full  range  of  pump  strengths  and  angular  tuning  (A:-axis). 
Figure  1  presents  a  succinct  overview  of  the  control  achieved  for  different  steering  angles  (proportional  to 
A;)  of  the  idealized  and  ramped  pump  systems  respectively.  For  these  plots,  the  external  pump  is  twice  the 
lasing  threshold  value,  the  feedback  strength  7  =  0.25,  and  the  feedback  delay  time  r  is  chosen  to  match 
the  selected  angle  (A;-value)  on  the  analytic  dispersion  curve.  Note  that  the  controlled  far-held  spectra 
are  extremely  sharp  in  k  and  oj  for  the  infinitely  wide  pump.  For  the  ramped  pump  case  in  (b),  there 
is  some  spectral  broadening  in  k  and  a;,  indicating  a  modulated  finite  support  traveling  wave  (Fig.  2.) 
rather  than  an  infinitely  extended  constant  amplitude  signal.  Besides,  the  controlled  state  is  no  longer  a 
traveling  wave  of  the  isolated  system  but  instead  is  a  solution  of  the  full  partial  differential  equation  with 
delay.  Its  wavenumber  (angle  tuning)  is  precise  but  the  frequency  is  off.  The  traveling  wave  emanates  on 
one  side  from  a  defect  (source)  and  is  absorbed  on  the  other  (a  sink),  as  shown  in  Fig.  2  for  r  =  2  and 
A;  =  4.  In  this  case  the  wave  travels  from  right  to  left  (since  the  corresponding  u  is  positive).  The  situation 
is  reversed  if  control  is  established  at  —k:  the  wave  then  travels  from  left  to  right.  Also  shown  on  this 
picture  is  beam  steering  from  A:  =  4  to  A:  =  7.  When  control  is  switched  from  one  wavenumber  to  the  other 
(around  t=0.27),  the  traveling  wave  at  A:  =  4  is  first  damped  out  since  it  is  no  longer  stabilized  by  the 
feedback  term.  Then  the  mode  at  A;  =  7  starts  growing  and  is  stabilized  around  t  =  6.8  (arbitrary  units). 
Our  observations  indicate  that  the  spatial  filter  is  the  critical  component  in  enabling  effective  control. 
The  time  delay,  while  less  critical,  proves  most  effective  when  matched  to  the  relevant  frequency  on  the 
nonlinear  dispersion  curve.  In  the  unramped  case,  we  have  seen  locking  to  an  intermediate  controlled  state 
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which  is  stable  but  not  the  desired  one  when  there  is  a  significant  frequency  mismatch  between  the  delay 
and  the  desired  frequency  on  the  dispersion  curve.  This  resembles  the  asymptotic  state  of  the  ramped 
system. 

The  control  scheme  has  been  extended  to  a  two-longitudinal  mode  model,  described  by  a  coupled  set  of 
CSH  equations[3].  This  study  shows  that  the  time  delay  (frequency  filter)  becomes  more  important  for 
successful  control  and  a  restriction  is  placed  on  the  feedback  control  strength. 
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Figure  1:  Left:  A  subset  of  controlled  k  —  u:  spectra  superimposed  on  the  analytic  dispersion  curve  (dashed 
line)  for  the  infinitely  extended  problem.  Right:  Similar  k  —  (jJ  spectra  for  the  ramped  pump. 
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Figure  2:  Left:  Modulated  traveling  waveform  and  envelope  for  the  controlled  finite  pump  case  with 
r  —  2,a  =  —5.  All  other  parameter  are  the  same.  Right:  Beam  steering  from  an  initial  controlled  signal 
at  A:  =  4  to  a  stable  signal  at  fc  =  7  at  the  same  parameter  values  as  in  the  left  figure. 
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Modeling  spatio-temporal  dynamics  of  broad  area/fiared  semiconductor 
laser  amplifiers 

J.K.  White,  J.G.  Mclnerney,  J.V.  Moloney,  S.W.  Koch 

Large  aperture  semiconductor  laser  amplifiers  are  extremely  promising  for  applications  combining  pulsed 
and  CW  power  levels  ~  1  -  10  W  with  high  spatial  coherence,  as  well  as  the  usual  advantages  of  semi¬ 
conductor  lasers.  This  project  relates  directly  to  modeling  ongoing  experiments  by  the  LIDA  group  at 
the  Phillips  Laboratory,  Kirtland  AFB,  Albuquerque.  The  model  incorporates  the  complete  spatiotempo- 
ral  dynamics  both  transverse  to  and  along  the  broad  area/flared  amplifier  or  laser  axis.  Existing  models 
extensively  used  in  the  general  literature  employ  the  Beam  Propagation  Method,  a  time  independent  ap¬ 
proach  which  in  the  specific  context  of  Icisers,  can  lead  to  ambiguity  in  the  interpretation  of  residts.  The 
present  approach  utilizes  a  look-up  table  for  the  gain  and  refractive  index  functions,  derived  from  the  basic 
microscopic  physics.  We  can  investigate  aU  intrinsically  time  dependent  behaviors  including  switch-on  of 
the  laser,  relaxation  oscillations,  uncontrollable  filamentation,  current  modulation,  injection-locking,  mir¬ 
ror  feedback  and  nonlinear  control,  arbitrary  current  injection  profiles  both  transverse  to  and  along  the 
laser  axis,  amplification  of  pulse  trains  in  a  broad  area/fiared  amplifier,  and  arbitrary  index/gain  guiding 
profiles  as  long  as  the  characteristic  time  scales  of  the  dynamics  exceeds  many  picoseconds.  Shorter  time 
scale  dynamics  requires  a  more  comprehensive  description  of  the  many-body  microscopic  physics  of  the 
interaction  of  light  with  semiconductor  materials.  In  addition,  some  issues  not  satisfactorily  dealt  with 
within  the  current  formulism  are  addressed  under  technical  report  Plasma  and  Lattice  Heating  Effects  in 
Quantum  Well  Lasers  and  VCSELs  and  in  [4]. 

The  basis  for  the  present  model  is  the  traveling  wave  amplifier  laser  rate  equations: 

ndE  dE  i  d'^E  ,  -r  \  c  Mn  r 

dN  ^  d^N  ,  J{x,z)  _  eocra5(iV,a)o)|pi2 

— Wv — 1^1 

where  the  symbols  have  their  usual  meanings.  Details  of  a  more  comprehensive  model  for  a  broad 
area/flared  laser  which  includes  counterpropagating  waves  and  boundary  conditions  have  already  been 
presented  in  [5] 

The  microscopic  physics,  including  many-body  effects  like  bandgap  renormalization  and  interband  carrier- 
carrier  scattering,  is  included  through  the  gain  function  g{N^oj^)  and  change  in  refractive  index 
both  being  determined  by  solving  the  semiconductor  Bloch  equations  in  the  quasi-equilibrium  approxima¬ 
tion  for  the  polarization  P: 
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The  Hopscotch  method  has  been  employed  to  solve  the  model  equations.  The  method  divides  the  com¬ 
putational  domain  into  a  checkerboard  of  alternating  even/odd  points.  At  each  time  step  the  odd/even 
points  are  advanced  explicitly  in  time  while  the  even/odd  points  are  then  advanced  implicitly.  The  next 
time  step  the  order  is  reversed.  This  method  has  the  advantage  of  properly  handling  the  beam  diffraction 
through  the  implicit  step  while  using  an  explicit  Lax-Wendrof  scheme  to  correctly  propagate  the  beam  in 
the  longitudinal  direction.  .Alternative  pseudo-spectral  methds  are  also  being  implemented  on  this  model. 
When  the  polarization  is  adiabatically  eliminated  from  the  Maxwell-Bloch  equations,  one  of  the  eigenvalues 
goes  marginally  unstable  for  large  values  of  the  transverse  wave  vector  k[l].  This  appears  as  a  large 
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oscillation  at  the  gridscale  of  the  computational  domain.  By  adding  a  small,  artificial  diffusion,  the  high-k 
oscillations  can  be  damped  out.  Care  must  be  given  that  the  diffusion  is  not  so  large  it  begins  to  dampen 
real  features  of  the  device.  A  more  physically  and  mathematically  consistent  approach  would  derive 
the  appropriate  order  parameter  equations  from  the  fuU  MajcweU- Semiconductor  Bloch  laser  equations,  a 
formidable  approach  currently  being  undertaken  as  part  of  the  present  contract. 

In  the  broad  area/flared  amplifier  case,  interaction  of  an  input  Gaussian  beam  with  the  sharp  edges  of  the 
current  stripe  leads  to  the  formation  of  “bat  ears”,  or  sharply  peaked  side  lobes.  This  implies  that  the 
sudden  drop  in  carrier  density  at  the  edge  of  the  current  stripe  is  responsible.  By  smoothing  this  edge  the 
“bat  ears”  can  be  reduced  (see  attached  figure).  A  linear  transverse  analysis  shows  that  a  perturbation 
related  to  the  first  and  second  derivatives  of  the  current  stripe  seeds  nonhnear  growth  of  the  “bat  ears” . 
Currently  we  are  investigating  multilongitudinal  modes  and  their  associated  dispersion  curves  in  broad 
area  semiconductor  lasers.  In  earlier  work  we  predicted  the  existence  of  parabolic  nonlinear  dispersion 
curves  (fc  vs  a;  )  for  a  free-running  broad  area  laser  undergoing  uncontrolled  filament  at  ion.  Such  dispersion 
curves  were  recenly  observed  experimentally  by  Dave  Bossert  at  the  Phillips  Lab/LIDA.  Our  current  goal 
is  to  correlate  our  theoretical  prediction  with  his  experimental  observation.  Eventually  we  hope  to  explore 
strategies  for  nonlinear -control  of  the  turbulent  output  of  these  devices.  Work  is  in  progress  on  a  joint 
publication  combining  theory  and  experiment. 

Modeling  the  spatiotemporal  dynamics  of  semiconductor  laser  amplifiers  presents  many  challenges.  By 
using  look-up  tables  to  incorporate  the  many-body  physics  and  artificial  diffusion  to  eliminate  the  high-k 
instability  we  have  created  a  model  that  captures  most  of  the  relevant  physics  in  these  devices.  Already 
we  have  shown  that  the  formation  of  “bat  ears”,  one  of  the  limiting  factors  of  power  in  the  central  lobe,  is 
a  result  of  the  input  pattern’s  interaction  with  the  edges  of  the  current  stripe 
Areas  of  future  study  include: 

•  FUamentation  dynamics 

•  Stability  under  pulsed  operation 

•  Derivation  of  complex  order  parameter  equations  for  the  semiconductor  problem  (see  J.  Lega,  J.V. 
Moloney,  and  A.C.  Newell,  “Universal  Description  of  Laser  Dynamics  Near  Threshold”,  Physica  D.) 

•  Developement  of  a  mean  field  code  to  be  run  in  “real  time”  on  a  Silicon  Graphics  “Power  Challenge” 
machine. 
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Full  Space-Time  Simulation  of  High  Brightness  Semiconductor  Lasers 
J.V.  Moloney,  C.Z.  Ning,  A.  Egan,  and  R.A.  Indik 

This  project  involves  direct  collaborative  links  to  Air  Force  Laboratory  scientists  (M.M.  Wright  and  D. 
Bossert  PL/LIDA)  and  to  industry  (R.  Patel,  S.  Srinivasan  Opto  Power  Corporation).  The  theoretical  basis 
for  this  model  involves  building  the  full  semiconductor  optical  response  function  from  a  microscopic  basis 
and  incorporating  this  in  the  laser  simulation  code.  Color  Plate  I  succinctly  summarizes  our  approach. 
Edge  emitting  wide  aperture  semiconductor  lasers  are  intrinsically  multi-moded  in  the  longitudinal  and 
transverse  dimensions.  Moreover,  the  carrier  density  and  field  intensities  display  large  nonuniformities 
both  along  and  transverse  to  the  device.  It  is  vital  therefore  to  take  account  of  the  gain  dispersion  (for 
mode  discrimination)  and  local  (in  space)  dependence  of  the  semiconductor  optical  response.  A  simulation 
model  which  takes  full  acount  of  the  above  features  has  recently  been  introduced  and  tested  against  the 
full  microscopic  many-body  model  for  pulse  propagation  in  a  bulk  GaAs  amplifier [1].  We  wUl  discuss 
recent  developements  in  the  microscopic  theory  which,  for  the  first  time,  allow  quantitative  comparison 
between  experimentally  measured  and  theoretically  computed  quantum  well  gain  spectra.  The  essence  of 
our  model,  discussed  in  detail  in  Ning  et  a/.[l],  is  that  the  local  gain  and  refractive  index  dispersion  as  a 
function  of  total  carrier  density  (N),  x(A,w),  extracted  from  the  microscopic  theory  or  experiment,  are 
fitted  with  multiple  Lorentzians; 
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where  Xo(A’)  is  a  “background”  contribution  which  is  frequency  independent.  Sq  =  Eg/h—uJc  with  Uc  being 
a  reference  frequency  and  Eg  the  bare  bandgap  energy.  The  number  of  Lorentzian  oscillators  required  for 
an  accurate  fit  wiU  depend  in  general  on  the  shape  of  the  gain/index  curves.  For  example,  quantum  weU 
structures  with  closely  spaced  subbands  may  require  more  terms  to  account  for  the  appearance  of  the 
higher  subband  peaks  in  the  gain  at  high  densities.  The  density  dependence  of  the  strength  A((A),  the 
width  Ti{N)  and  the  position  Si(N)  in  the  frequency  domain  of  the  individual  Lorentzian  oscillators  allows 
the  laser  to  respond  on  the  fly  to  local  changes  in  N  along  and  across  the  structure  as  a  function  of  time. 
It  is  important  to  stress  that  both  semiconductor  material  dimensionality  effects  (3D,  2D),  strain  etc  and 
the  full  many-body  effects  such  as  bandgap  renormalization.  Coulomb  screening,  etc.  are  all  included  in 
the  microscopically  computed  gain  and  refractive  index  spectra  x(A,u)). 

From  each  oscillator  can  construct  a  corresponding  polarization  Pj{u))  through  the  relation 
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eoe,Aj{N)E(io) 

Tj{N)  +  i(So+u;-6g{Ny 


(j  =  0,1,2...) 


(2) 


where  quantities  with  a  "represent  their  Fourier  transforms.  E  is  the  laser  field  amplitude.  Inverting  the 
Fourier  transform  for  (j  =  1,2, ...),  we  obtain 


=  {-TjiN)  +  i[6o  -  6jiN)]}Pjit)  -  i€oe,MN)E{t),  (3) 

while  for  (j  =  0)  we  simply  have 


Po(iV)  =  eoe,Xo(N)Eit),  (4) 

because  no  frequency  dependence  of  Xo(A)  is  assumed.  The  assumption  that  the  carrier  density  N  is 
frozen  on  the  time  scale  of  the  polarization  dynamics  is  a  good  approximation  for  a  laser. 

Both  the  forward  and  backward  field  and  carrier  density  equations  must  now  be  appropriately  modified  to 
couple  to  the  locally  linear  oscillators. 
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where  Ug  is  the  group  index,  and  K  =  (^c'^bl^  is  the  wavevector  in  the  medium  of  background  index  tij. 
The  ?/,  J,  e,  w  are  respectively,  the  quantum  efficiency,  pumping  current,  electron  charge,  and  active  region 
width.  Once  a  laser  structure  is  specified,  the  corresponding  gain  and  index  spectra  can  be  computed  and 
incorporated  into  the  counterpropagating  field  code. 

As  a  test  of  this  model  we  have  carried  out  a  full  space-time  simulation  of  the  field  and  carrier  density 
build-up  from  noise  in  a  MOPA  device.  The  device  structure  consists  of  two  sections,  a  DBR  laser  Master 
Oscillator  (MO)  of  total  length  1.575  mm  and  a  linearly  expanding  Power  Amplifier  (PA)  of  length  2  mm. 
The  MO  has  a  750  fim  pumped  gain  section,  a  750  fim  DBR  end  distributed  reflector  and  a  shorter  75  fim 
DBR  reflector  at  the  PA  end.  The  PA  has  a  separate  electrical  contact.  The  passive  semiconductor  under 
both  DBR  sections  needs  to  be  at  least,  partially  bleached,  before  the  forward  and  backward  waves  can 
see  the  gratings.  This  is  reflected  in  the  rather  sharp  turn-on  of  the  device  as  the  MO  current  is  increased. 
Another  important  issue  is  the  weak  reflectivity  from  the  PA  output  facet.  This  is  known  to  be  detrimental 
to  the  stable  operation  of  the  device.  We  have  studied  the  spatiotemporal  dynamics  of  the  laser  turn-on 
for  a  PA  facet  output  reflectivity  Rp  =  0%,  0.05%  and  0.1%. 

Surface  plots  of  the  carrier  density  and  field  intensity  are  shown  above  and  below  the  device  geometry  in 
Color  Plate  II,  for  stable  operation  when  when  Rf  =  0.05%.  Here  the  MO  current  is  Imo  =  100  mA  and 
the  PA  current  is  Ip  a  =  1.5  A.  The  forward  field  at  the  beginning  of  the  PA  does  not  follow  the  diflFraction 
pattern  predicted  for  a  finite  aperture  of  4  //m  as  strong  gain  and  index  guiding  provide  a  nonlinear 
modification.  The  lack  of  optical  field  between  the  intensely  illuminated  region  and  the  straight  edge  of  the 
PA  leads  to  a  pronounced  concentration  of  unsaturated  carriers  in  this  region.  These  unsaturated  carriers 
wiU  be  available  for  gain  if  there  is  any  finite  amplified  back  reflection  in  the  PA.  We  also  note  that  the 
MO  output  intensity  is  insufficient  to  saturate  the  carriers  just  beyond  the  entrance  to  the  PA. 

Turning  on  a  finite  facet  reflectivity,  one  immediately  sees  evidence  for  an  external  resonator  behavior  [2]. 
The  important  new  feature  is  the  finite  back-reflected  field  which  becomes  strongly  amplified  in  the  region 
of  unsaturated  carriers  and  leaks  through  the  lossy  75  /xm  DBR  into  the  MO.  The  backward  field  also 
displays  some  fine  features  in  the  section  of  MO  close  to  the  interface  with  the  PA,  which  is  a  signature 
of  the  multimode  behavior.  As  a  result,  the  MOPA  becomes  unstable.  The  single  mode  integrity  of 
the  MO  is  now  destroyed  and  one  can  observe  multi-longitudinal  mode  oscillations  associated  with  this 
section  at  frequencies  typically  in  the  range  of  30  —  40  GHz.  We  have  also  observed  multi-longitudinal 
mode  beats  at  around  16  GHz  corresponding  to  the  PA  section.  The  transverse  profile  remains  essentially 
unchanged  in  shape  in  this  unstable  regime.  At  much  higher  PA  pumping  currents  we  have  observed 
dynamic  filamentation  instabilities  in  addition  to  the  multi-longitudinal  mode  beatings. 

Preliminary  studies  of  current  modulation  of  the  MOPA  device  have  already  been  carried  out.  Modulation 
of  MOPAs  is  an  important  consideration  in  potential  communications  applications.  Experimental  investi¬ 
gations  of  integrated  MOPAs  under  current  modulation  (CM)  and  electro-absorption  modulation  (EAM) 
of  the  MO  have  taken  place,  and  have  demonstrated  some  of  the  advantages  and  limitations  that  may 
exist  in  both  schemes.  For  CM  the  main  limitation  to  achieving  high  power  is  the  tendency  for  the  PA 
to  lase  independently  of  the  MO,  when  the  MO  current  remains  at  or  below  its  threshold  value  for  longer 
than  the  turn-on  time  delay  in  the  PA.  Another  problem  is  that  transverse  mode  instabilities  may  arise 
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at  liigh  PA  currents.  For  example,  Figure  1(a)  shows  the  temporal  evolution  of  the  output  power  from 
the  MOPA  under  sinusoidal  current  modulation  of  the  MO  at  4  GHz,  close  to  the  threshold  for  transverse 
mode  instabilites  in  the  device.  Under  DC  bias  conditions  the  output  power  is  stable  and  constant,  but  it 
becomes  irregular  and  unstable  under  modulation.  Our  model  permits  us  to  resolve  the  transverse  mode 
behavior  (Figure  1(c))  and  determine  how  the  profile  changes  in  the  course  of  the  emission  of  a  pulse. 

We  can  also  explore  the  spatio-temporal  dynamics  of  the  MOPA  under  EAM.  An  extra  section  is  added  to 
the  MOPA  between  the  MO  and  PA,  which  can  be  switched  between  a  nearly  transparent  “on”  state  and  a 
nearly  opaque  “off”  state  with  a  reverse  bias.  This  arrangement  suppresses  the  tendency  for  the  PA  to  lase 
independently  of  the  MO.  However,  during  the  “off”  states,  the  carrier  density  in  the  PA  rises  to  a  high 
level,  so  that  when  the  EAM  is  switched  to  the  “on”  state  considerable  feedback  occurs  and  undamping 
of  the  relaxation  oscillation  in  the  MO  occurs.  This  places  a  limitation  on  the  PA  operating  current  and 
hence  on  the  output  power. 

Our  simulation  results  suggest  a  modification  in  the  design  which  should  increase  this  window  of  stable 
operation,  namely  designing  the  PA  current  pump  as  a  trumpet  flare  matching  the  actual  beam  divergence 
from  the  MO  output  facet.  Color  Plate  III  contrasts  the  new  design  with  the  existing  linear  flared  device. 
Discussions  are  currently  underway  with  Opto  Power  Corporation  to  build  a  prototype  high-brightness 
source  based  on  this  design.  We  emphasize  that  the  model  used  here  and  developed  in  detail  in  reference 
[1],  is  not  confined  to  the  MOPA  structure  but  can  be  used  to  study  single  and  multi- section  lasers  of 
arbitrary  geometry.  We  are  in  the  process  of  including  thermal  processes  at  the  level  of  ultrafast  plasma 
and  lattice  heating  through  to  quasistatic  bulk  heating  effects. 
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Figure  Captions 

Figure  1 :  For  IpA  =  2.75A,  (a)  the  total  power  output  from  the  device  and  the  power  contained  to  the  Full 
Width  at  Half  Maximum  (FWHM),  (b)  the  modulated  current  applied  to  the  MO  and  (c)  the  evolution 
of  the  FWHM  of  the  central  lobe  of  the  farfield  intensity  profile. 
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Figure  1 :  For  Ip  a  -  2.75A,  (a)  the  total  power  output  from  the  device  and  the  power  contained  to  the  M 
Width  at  Half  Maximum  (FWHM),  (b)  the  modulated  current  applied  to  the  MO  and  (c)  the  evolution 
of  the  FWHM  of  the  central  lobe  of  the  farfield  intensity  profile. 
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Plasma  and  Lattice  Heating  Effects  in  Quantum  Well  Lasers  and  VCSELs 
C.Z,  Ning,  R.  Indik,  J.V.  Moloney 

Summary 

Heating  or  thermal  effects  ctre  a  serious  problem  for  semiconductor  lasers  in  general  and  for  some  state-of- 
the-art  devices,  such  as  the  Vertical- Cavity  Surface- Emitting  Lasers  (VCSELs)  in  particular.  There  is  much 
experimental  evidence  showing  that  heating  can  influence  the  device  operation  significantly.  Furthermore, 
both  theoretical  and  expeiimental  investigations  indicate  that  the  lattice  and  plasma  can  each  attain 
individual  temperatures  that  are  generally  different  from  the  ambient  temperature  during  laser  operation. 
It  is  thus  of  crucial  importance  to  develop  a  self-consistent  theory  of  heating  in  semiconductor  lasers 
including  lattice  and  plasma  temperatures  as  independent  dynamical  variables,  in  order  to  understand 
the  underlying  physics  and  to  help  design  semiconductor  laser  devices.  We  have  developed  such  a  self- 
consistent  theory  and  applied  it  to  CW  operation  of  VCSELs.  Many  phenomena  observed  experimentally 
can  be  explained  using  our  approach.  These  include  the  increase  of  CW-density  and  change  of  frequency 
shift  with  pumping,  and  thermal  switch  off.  Furthermore  our  approach  predicts  a  novel  threshold  bistability 
of  VCSELs.  The  well-known  gain  nonlinearity  can  also  be  justified  or  improved  using  our  present  theory. 


The  model 

Our  starting  point  is  the  free-carrier  model  for  the  semiconductor  laser  with  k-resolved  polarizations  adi- 
abaticaUy  eliminated.  This  set  of  equations  is  generally  equivalent  to  an  infinite  hierarchy  of  moment 
equations.  We  cut  off  the  infinite  hierarchy  at  the  second  order  moment  and  close  the  hierarchy  by  as¬ 
suming  that  the  k-resolved  densities  obey  the  Fermi  distribution.  This  leads  to  two  equations  for  the  total 
carrier  density  (N)  and  energy  (W).  The  energy  equation  is  transformed  to  an  equation  for  plasma  tem¬ 
perature  {Tp).  The  lattice  temperature  equation  is  established  phenomenologically  by  considering  different 
cooling  and  heating  mechanisms  of  the  lattice.  The  coupling  among  different  parts  of  the  total  system  is 
schematically  shown  in  Fig.l.  The  whole  set  of  equations  are  given  as  follows  (with  A  being  the  complex 
amplitude  of  the  field): 
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Plasma  heating 

Plasma  heating  in  CW  operation  is  studied  by  fixing  the  lattice  temperature  at  300  K  and  solving  the  set 
of  nonlinear  equations  (1-3)  self-consistently.  We  observe  significant  plasma  heating  effects  including  the 
dependence  of  the  cw-carrier  density  and  frequency  shift  on  pumping.  An  example  of  the  CW  solutions  is 
shown  in  Fig.2. 


Gain  nonlinearity 

Gain  nonlinearity  (or  gain  compression)  means  a  nonlinear  dependence  of  the  gain  on  the  laser  intensity  in 
the  semiconductor  rate  equations.  This  nonlinear  gain  is  very  important  in  reconciling  the  disagreement 
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between  experiments  on  diiectly  modulated  lasers  and  the  simple  rate  equation  (without  nonlinear  gain) 
predictions.  Here  we  show  that  plasma  heating  leads  to  a  nonlinear  gain:  Starting  from  equations  (1-3), 
we  eliminate  the  plasma  temperature  adiabatically.  This  leads  to  an  elevated  plasma  temperature  relative 
to  the  lattice  temperature.  Formally  we  have: 

tp  =  0  T,  =  Tj,{N,\A?)  (5) 

Substitute  this  stationary  plasma  temperature  in  to  equations  (1)  and  (2)  we  have: 


X{T,,  N)-^x  {Tp{N,  \A\%  N)  =  x(^,  |  Ap)  (6) 

By  solving  the  plasma  temperature  equation  we  show  that  the  resulting  nonlinear  gain  can  be  fitted  by 
the  following  non-perturbative  form 


-Imx{NM\") 


(7) 


The  fitting  to  (7)  (diamonds)  and  the  solution  to  (5)  (solid  line)  for  different  heat  transfer  constants  are 
shown  in  Fig.3 


Thermal  switch-off 

In  experiments  one  often  observes  a  strong  saturation  or  even  switch-off  of  the  laser  output  at  quite  mod¬ 
erate  pumping  strength.  To  understand  this  phenomenon  theoretically,  we  solve  numerically  the  whole  set 
of  equations  (1-4).  Examples  from  our  CW-solutions  are  shown  in  Fig.4 


Thershold  bistability 

Under  appropriate  conditions  the  threshold  of  VCSELs  can  become  bistable  as  shown  in  Fig.5.  This 
rather  unexpected  phenomenon  results  from  the  complicated  interplay  between  lasing  on  the  one  hand 
and  plasma  and  lattice  heating  on  the  other.  The  occurrence  of  this  phenomenon  depends  on  the  ambient 
temperature  and  the  rate  of  heat  transfer  to  the  device  substrate.  Currently  we  are  collaborating  with 
a  group  of  experimentalists  at  University  of  New  Mexico  to  compare  quantitatively  this  predication  with 
experiments. 

Future  directions 

The  research  reported  above  needs  to  be  extended  in  several  directions  and  some  related  work  stiU  needs 
to  be  done.  The  following  directions  are  of  special  interest: 

•  Inclusion  of  many-body  interactions  in  our  present  formalism.  This  would  be  very  important  if  a 
quantitative  comparison  with  experiments  is  intended,  because  of  the  band  gap  renormalization  due 
to  these  interactions.  This  research  is  currently  being  carried  out  in  a  collaboration  with  Stephan 
Koch  of  the  University  of  Marburg  in  Germany. 

•  Investigations  of  the  transverse  structure  in  VCSELs,  especially  under  the  transverse  heat  conduction 
and  carrier  diffusion. 

•  Dynamics  and  patterns  in  coupled  VCSEL-arrays  using  a  parameterized  coupled  rate-equation  ap¬ 
proach. 
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•  The  influence  of  plasria  and  lattice  heating  on  modulation  bandwidth. 


•  A  noise  theory  including  thermal  effects  that  could  provide  us  with  a  correct  linewidth  enhancement 
factor. 
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Polarization  Patterns  in  Wide  Aperture  3-level  and  Semiconductor  Lasers 

Polarization  pattern  forming  instabilities  in  lasers 
Q.  Feng,  M.  San  Miguel.,  J.V.  Moloney  and  A.C.  Newell 

Transverse  pattern  formation  in  lasers  has  attracted  much  attention  in  recent  years  and  a  rich  variety 
of  pattern  forming  instabilities  has  been  found.  Most  of  studies  have  concentrated  on  systems  in  which 
the  direction  of  the  electric  field  in  the  laser  cavity  is  fixed  (by,  e.g.  Brewster  windows).  In  this  study 
we  consider  the  case  where  the  degree  of  freedom  of  the  electric  field  polarization  is  unconstrained.  We 
analyze  the  stability  of  solutions  corresponding  to  transverse  patterns:  linearly  polarized  traveling  waves 
and  polarized  standing  waves.  The  interplay  between  the  polarization  and  transverse  effects  leads  to  new 
instabilities  whose  natures  are  identified  with  amplitide  equations.  Here  we  report  results  for  the  positive 
detuning  case  (results  for  negative  detuning  have  been  reported  elsewhere[l]). 

The  system  we  considered  is  a  wide  aperture,  single  longitudinal  mode,  ring  cavity  laser  with  transverse  flat 
end  reflectors.  The  basic  equations  are  the  vectorial  MaweU- Bloch  equationsgoverning  space-time  evolution 
of  the  slowly  varying  complex  envelopes  of  the  electric  field  in  the  laser  cavity,  E±  =  ^{Ex  ±  iEy),  which 
are  the  left-  and  right-circularly  polarized  components.  For  positive  detuning,  lasing  modes  with  finite 
transverse  wavenumbers  lose  stability  as  the  pump  (r)  is  increased  above  threshold  (rc).  From  the  vectorial 
MaxweU-Bloch  equations  we  derived  the  amplitude  equations  which  are  valid  near  threshold: 
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where  Zi  (2^2,  and  z^)  is  the  amplitude  of  right-  (right-,  left-  and  left-)  circularly  polarized,  traveling 
wave  traveling  in  +x  {—x,  —x  and  -fa:)  direction,  respectively.  The  coefficients  a  and  /?  are  related  to  the 
material  parameters:  a  =  6“^,  (3  =  -f  c”^)/2,  where  b  is  the  decay  rate  of  the  population  inversion 
between  the  atomic  levels  (J  =  1,  =  ±1)  and  (J  =  0),  c  the  decay  rate  of  the  coherence  between  the 

levels  (J  =  1,  =  ±1)  (they  are  measured  in  units  of  material  polarization  decay  rate).  Usually  c>  b 

and  we  wiU  restrict  ourself  to  this  case  through  this  study. 

The  Amplitude  equations  admit  linearly  polarized  traveling  wave  solutions, 

{zi,Z2,Z3,Z4)  =  (t;,0,0,n)e‘(‘^®"‘'*\ 


A  linear  stability  analysis  determines  the  stability  wavenumber  \Q\  band  (for  given  e  =  r  —  r^).  We  found 
four  types  of  instabilities,  each  of  which  sets  in  at  wavenumber  \Q\  >  \Qi\^  i  =  1, 2, 3, 4,  respectively,  where 


ijL  —  1_L 

3^2’  ^2  5  +  7c^2’  ^3  2^2’ 


Ql  = 


c  —  b  € 
4c  ’ 


The  instability  occuring  at  Qi  is  the  well-known  Eckhaus  instability  (the  typical  --rule),  which  is  caused 
here  by  the  long- wavelength  perturbation  in  the  sum  of  the  phases  of  Zi  and  2:4,  the  right-  and  left- circularly 
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Figure  3:  Stability  diagrams  for  positive  detuning.  The  open  circles,  filled  circles,  crosses  and  asterisks 
correspond  to  Qi,  Q2,  Qa  ^^d  Qs,  respectively  (see  text).  The  dotted  line  is  the  neutral  curve.  The  solid 
and  dashed  lines  are  stability  boundaries  obtained  from  direct  stability  analysis  based  on  the  vectorial 
MaxwelTBloch  equations. 

polarized  traveling  waves  (traveling  in  the  same  direction).  The  Q2-instability  is  the  one  due  to  the  long- 
wavelength  perturbation  in  the  phase  difference  of  Zi  and  Z4.  The  Qs-instability  is  the  amplitude  instability 
due  the  interaction  with  a  wave,  with  the  same  polarization,  traveling  in  the  opposite  direction[2].  The 
Q 4-instability  is  another  amplitude  instability  which  is  caused  by  an  orthogonally  polarized,  oppositely 
traveling  wave.  Both  Q2  and  Q4  depend  on  relative  magnitudes  of  the  decay  rates:  Q2  0*^  and  Q4  O'*" 
as  c  This  means  that  the  stable  band  disappears  as  c  — ^  b'^.  Near  threshold  the  -instability 

limits  the  stability  band,  while  far  above  threshold  it  is  the  amplitude  instability  originating  from  the 
Q4-instability  that  limits  the  stability  band  (see  Fig.l).  The  instabilities  in  the  far  above  threshold  region 
are  obtained  from  the  numerical  stability  analysis  based  on  the  vectorial  MaxweU-Bloch  equations. 

The  amplitude  equations  have  another  type  of  stable  solutions,  the  polarized  standing  waves, 

{zuZ2,Z3,Z4)  = 

which  correspond  to  standing  waves  in  both  x  and  y  components  of  the  electric  field.  We  have  also 
studied  the  stability  of  these  solutions  and  found  four  types  of  instabilities:  two  phase  instabilities  and  two 
amplitude  instabilities.  Among  the  two  phase  instabilities,  however,  the  Eckhaus  |-rule-instability  is  now 
caused  by  the  long- wavelength  perturbation  in  the  phase  difference  of  Zi  and  Zs,  while  the  instability  which 
sets  in  at  Q 2  is  caused  by  the  long- wavelength  perturbation  in  the  sum  of  the  phases.  We  have  performed  a 
series  of  numerical  simulation  on  the  vectorial  MaxweU-Bloch  equations  with  different  transverse  boundary 
conditions.  The  linearly  polarized  traveling  waves  tend  to  be  destabilized,  while  the  polarized  standing 
waves  may  be  more  robust  to  survive  under  reaUstic  boundaries. 
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Light  polarization  dynamics  in  surface  emitting  semiconductor  lasers 
M.  San  Miguel,  Q.  Feng  and  J.V.  Moloney 

The  basic  modeling  of  semiconductor  laser  dynamics  is  provided  by  rate  equations  and  their  generalization 
to  include  the  dynamics  oi'  the  phase  of  the  electric  field.  The  rate  equations  take  for  granted  a  fixed 
direction  of  polarization.  Although  they  have  been  very  useful  both  in  basic  and  applied  research  in 
semiconductor  laser  devices,  no  such  laser  modeling  seems  to  be  available  to  describe  the  polarization  of 
the  electric  field.  On  the  other  hand,  a  number  of  polarization-sensitive  applications  of  semiconductor 
laser  devices  require  detailed  polarization  control,  and  polarization  coexistence  between  orthogonal  modes 
(polarization  partition)  and  polarization  switching  have  been  observed.  The  purpose  of  this  work[l]  is  to 
introduce  and  analyze  a  basic  model  of  semiconductor  dynamics  in  which  the  key  aspects  of  polarization 
dynamics  of  a  semiconductor  laser  can  be  studied.  The  semiconductor  engineering  involving  polarization 
control  would  benefit  from  a  basic  modeling  of  the  key  physical  issues. 

We  introduce  a  four-level  model  which  takes  into  account  the  spin  sublevels  of  the  conduction  and  valence 
band  and  therefore  allows  us  to  consider  different  polarizations  of  light  associated  with  transitions  between 
different  spin  sublevels  (see  the  included  figure  for  a  schematic  illustration).  Generalized  rate  equations 
which  include  polarization  degrees  of  freedom  are  obtained  from  the  four-level  model.  Analysis  of  these 
equations  privodes  physical  important  insight  into  the  nature  of  polarization  instabilities  in  surface  emitting 
semiconductor  lasers.  In  particular,  we  show  that  the  type  dynamical  response  of  the  polarization  degrees 
of  freedom  is  linked  to  the  relative  time  scale  of  spontaneous-emission  and  spin-relaxation  processes. 

The  polarization  dynamics  is  closely  interrelated  with  transverse  effects  in  lasers  with  large  Fresnel  number. 
In  order  to  study  general  issues  of  the  coupling  of  the  polarization  and  the  transverse  degrees  of  freedom, 
we  have  included  in  our  vectorial  rate  equations  a  general  modeling  of  transverse  effects  for  broad  area 
lasers.  Near  threshold  we  derive  an  amplitude  equation  description  which  avoids  the  unphysical  high- 
wavenumber  instability  resulted  from  the  adiabatic  elimination  of  the  dipole  polarization  variables  in 
the  usually  approach[2,  3].  We  predict  the  existence  of  stable  transverse  spatially  homogeneous  intensity 
outputs,  corresponding  to  off-axis  emissions,  with  arbitrary  direction  of  linear  polarization  in  the  transverse 
plane.  The  stability  of  the  off-axis  emission  solutions  to  long-wavelength  perturbations  is  investigated  and, 
in  addition  to  the  Eckhaus  instability  associated  with  a  global  phase,  we  predict  a  new  type  of  polarization 
instability  associated  with  a  relative  pha^e  of  the  complex  field  vector.  We  also  explore  the  role  of  laser 
cavity  anisotropies  and  their  interplay  with  transverse  effects.  Linear  phase  anisotropies  associated  with 
birefrigence  stabilize  two  prefered  orthogonal  directions  for  off-axis  linearly  polarized  light.  An  important 
conclusion  is  that  transverse  effects  discriminate  between  these  two  prefered  directions  so  that,  for  given 
sign  of  the  anisotropy,  they  extend  the  range  of  stability  of  the  x-polarized  emission  while  the  stability 
range  of  the  ^/-polarized  solution  shrinks. 
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1D-3D  Vector  Maxwell  Solvers 
R.G.  Flesch  and  J.V.  Moloney 

We  have  developed  a  suite  of  modular  one,  two  and  three-dimensional  vector  Maxwell  solvers  to  be  used 
in  the  study  of  the  propagciion  of  ultrashort  (<  100  fs)  pulses  in  a  variety  of  nonlinear  materials.  These 
codes  serve  two  purposes: 

(1)  To  solve  problems  where  a  vector  Maxwell  description  is  essential. 

(2)  To  validate  envelope  theories  (and  their  corrections)  is  regimes  where  the  envelope  assumption  becomes 
suspect. 

The  former  category  includes  ultrashort  intense  pulses  propagating  through  structures  with  wavelength 
and  sub-wavelength  dimensional  inhomogeneities.  Examples  are:  periodic  linear /nonlinear  index  gratings, 
Quantum-Well  stacks,  quantum  wires,  semiconductor  micro-cavities,  linear/nonlinear  surfaces  and  the 
vicinity  of  the  critical  focus  of  a  collapsing  beam.  In  ID  many  of  the  simulations  can  be  carried  out 
in  real  time  on  our  new  Power  Challenge  L  machine.  This  offers  the  opportunity  to  develop  a  CAD- 
like  interactive  computing  and  graphical  visualization  environment  which  will  significantly  enhance  our 
efficiency  in  modeling  a  wide  variety  of  nonlinear  optical  structures. 

Maxwell’s  equations  have  the  following  form: 


dt 


-V  X  E 


dt 


— V  X  £ 
Mo 


(1) 


The  equations  ioi  E  and  B  are  closed  by  adding  the  constitutive  relation 


D  —  eoSooE  +  Pl  +  Pnl  (2) 

which  allows  for  both  a  linear  and  nonlinear  polarizations  which  we  may  write  in  general  as 

Pl  =  Ol{D,E)  (3) 

Pnl  =  •  (4) 

The  complete  solution  of  the  coupled  system  consists  of  the  following  three  steps: 

1.  Given  values  for  E  and  P,  obtain  updates  for  B  and  D.  The  Maxwell  equations  are  solved  by  the 
FT-TD  method. 

2.  Using  the  values  from  step  1,  calculate  the  linear  and  nonlinear  polarizations. 

3.  Obtain  a  new  value  for  E  from  the  constitutive  relationship  given  jD,  Pl  and  Pjvn  from  steps  1  and 

2. 

To  date  we  have  used  a  linear  Lorentz  medium  for  the  linear  polarization  vector  and  a  mixed  Kerr  and 
Raman  response  for  the  nonlinear  polarizations.  To  verify  and  debug  the  code  we  have  investigated  the 
following  situations: 

•  Reflection  of  a  pulse  from  a  layered  set  of  dielectrics. 

•  Soliton  formation  for  very  short  (<  100  fs  FWHM)  pulses. 

•  Propagation  of  ‘‘superluminal”  pulses, 

•  Precursor  phenomena. 

•  Carrier  wave  shocks. 
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•  Paraxial  corrections  to  the  linear  focusing  (in  vacuum)  of  a  very  narrow  (two  wavelengths  FWHM) 
Gaussian  pulse. 

•  Nonlinear  self  focusing  of  a  Gaussian  pulse. 


The  shocking  of  the  carrier  wave  is  of  particular  interest  and  to  our  knowledge  has  not  yet  been  addressed 
in  the  literature.  Carrier  shocking,  as  depicted  in  Fig  1,  is  due  to  the  nonlinear  dependence  of  the  wave 
velocity  on  the  amplitude.  Using  a  simple  analysis  for  a  medium  without  dispersion  one  can  predict  a 
shock  time  for  the  carrier  which  is  independent  of  the  pulse  width  for  pulses  which  contain  at  least  two 
wavelengths.  Assuming  a  Kerr  nonlinearity  Pjvr  =  •  E)E  the  estimated  time  to  breaking  of  the 

carrier  wave  is  given  by 


Tb  =  - 


2  T 


opt 


(1  +  A.hx^^'^Elf'^ 


3  X^^^El  ^(1  +  3x(^)Po')(l  +  ^X^^'^El)  ' 


(5) 


In  Fig.  2  we  present  results  of  the  numerically-determined  breaking  times  which  agree  rather  well  with  the 
results  in  Eq.  5.  Included  in  this  figure  is  the  curve  (dots)  of  envelope  breaking  times  for  an  optical  pulse 
of  300  fs  duration. 

To  model  a  more  realistic  situation  we  add  a  small  amount  of  dispersion.  This  introduces  a  mismatch 
between  the  phase  and  group  velocities  of  both  the  fundamental  and  harmonic  freqencies  which  in  most 
cases  prevents  carrier  wave  shock  formation.  As  a  pulse  enters  the  nonlinear  medium  we  see  the  carrier 
steepen,  followed  by  the  harmonic  components  separating  according  to  their  phase  velocities.  The  presence 
of  dispersion  and  hence  nonzero  GVD  allows  for  envelope  shock  formation  which  has  also  been  observed. 
Carrier  shock  formation  is  of  fundamental  importance  in  the  study  of  the  critical  collapse  of  a  Gaussian 
packet  in  three  spatial  dimensions  because  the  high  field  amplitudes  which  occur  during  the  collapse 
can  give  rise  to  carrier  shock.  The  dynamics  is  extremely  rich  due  to  the  complex  interplay  of  the  carrier 
shocks,  envelope  shocks  and  self-focusing.  For  a  full  3D  simulation  of  the  coupled  system  we  require  several 
Gigabytes  of  memory  and  the  equivalent  of  hours  of  Cray  CPU  time.  The  large  memory  requirements  are 
due  to  fact  that  the  carrier  frequency  (and  possibly  third  harmonics)  of  the  pulse  must  be  well  enough 
resolved  to  avoid  numerical  dispersion  from  overwhelming  any  actual  dispersion  in  the  problem.  The 
code  has  been  implemented  and  optimized  for  the  Vicksburg  C90  and  the  AHPCRC  CM5.  We  anticipate 
developing  a  PVM  version  which  wiU  run  on  the  Mauii  SP2 

As  an  example  of  our  three-dimensional  results,  we  present  preliminary  results  of  the  linear  and  nonlinear 
focusing  of  a  Gaussian  beam  in  a  medium  without  dispersion.  In  the  linear  focusing  case  we  see  in  Fig.  3 
that  in  addition  to  the  original  x  polarization  (a),  a  longitudinal  (b)  and  additional  transverse  component 
(c)  develop  as  the  pulse  approaches  the  linear  focal  point.  The  transverse  structure  of  each  of  these 
polarizations  is  shown  in  Fig.  4.  These  polarizations  represent  corrections  to  the  paraxial  approximation 
in  which  only  the  initial  transverse  polarization  is  of  order  one.  To  study  critical  collapse  we  wish  to 
concentrate  on  the  nonlinear  effects  and  hence  we  choose  a  beam  which  is  much  wider  thus  initially  avoid 
strong  paraxial  corrections.  As  the  pulse  propagates  further  the  width  becomes  too  narrow  and  our 
transverse  resolution  is  too  small.  As  in  the  one  dimensional  dispersionless  case  we  also  observe  carrier 
shock  here.  Although  the  addition  of  some  dispersion  prevents  the  carrier  shock  from  forming  and  allows 
us  to  follow  the  focusing  further,  we  have  not  yet  been  able  to  observe  the  inhibition  of  self-focusing  due 
to  the  GVD. 
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Figure  3 


Dynamic  Nonlinear  Optical  Skin  Effect 
W.  Forysiak,  E.M.  Wright  and  J.  V.  Moloney 

The  propagation  of  an  intense  optical  field  incident  from  air  on  a  saturable  absorber  can  give  rise  to  a 
self-reflected  wave,  if  the  absorption  is  sufficiently  large.  The  reflected  wave  arises  at  a  wavelength-scale, 
spatial  transition  between  the  saturated  and  unsaturated  regions  of  the  absorber,  akin  to  the  skin  effect 
at  a  linear  boundary,  and  is  self-reflected  in  the  sense  that  the  transition  region  is  induced  by  the  incident 
field  itself. 

The  FD-TD  method  discretises  the  differential  form  of  Maxwell’s  equations  directly  and  allows  one  to 
determine  the  evolution  of  an  optical  field  in  a  nonlinear  medium,  subject  to  the  given  constitutive  relations 
between  the  electric  field  and  polarisation,  and  without  recourse  to  the  SVEA.  We  examine  pulse  shaping 
in  the  near-field,  close  to  the  absorber  boundary,  and  predict  a  new  nonlinear  focusing  effect,  which  we 
attribute  to  the  formation  of  a  transient,  focusing  mirror  in  the  absorber.  As  the  incident  pulse  impinges 
on  and  strongly  saturates  the  absorber,  it  excites  a  moving  reflection  front  [1]  which  is  shaped  according  to 
the  transverse  profile  of  the  incident  pulse.  If  the  incident  pulse  transverse  profile  is  bell-shaped,  so  too  is 
the  resulting  mirror,  and  the  reflected  pulse  is  focused  according  to  the  waist  and  intensity  of  the  incident 
pulse.  In  addition  to  being  re-shaped,  the  reflected  pulse  is  spectrally  broadened  and  red-shifted  due  to 
the  Doppler  effect  at  the  moving  mirror  [1]. 

We  consider  the  time-dependent  propagation  of  a  2-D  transverse  electric  (TE)  polarised  pulse,  in  which  the 
electric  field  is  polarized  along  the  y-axis  and  also  assumed  uniform  along  that  axis,  E{r,t)  =  yEy{x,z,t). 
Then  the  Maxwell’s  equations  for  the  electric  and  magnetic  field  quantities,  Ey,  and  are, 

dP^  _  1  fdB^  dBA 

dt  dz  '  dt  dx  ’  dt  piQ\dz  dx  )  ' 

where  the  z-axis  is  the  propagation  direction,  and  the  x-axis  is  the  transverse  direction.  The  nonlinear 
optical  response  of  the  saturable  absorber  is  included  using  a  two-level  model  via  the  constitutive  relation, 
Dy  =  e^Ey  -h  Py,  whepo  the  macroscopic  polarisation,  Py  =  Np(p2i  -f  c.c.),  is  determined  by  the  Bloch 
equations, 

+  i'fi  +  ioJ2i)p2i  =  i^n,  ^  +  jiin  -  1)  =  2i^{p2i  -  p*2i).  (2) 

Here,  N  is  the  density,  p2i  is  the  off-diagonal  density  matrix  element,  n  =  {pn  —  P22)  is  the  population 
difference,  U21  is  the  transition  frequency,  p  is  the  dipole  moment,  and  71  —  l/^i  and  72  =  1/^2  ^-re  the 
population  and  polarization  damping  constants. 

The  conditions  for  self-reflection  of  a  continuous  plane-wave  incident  field  require  that  the  normalised 
parameters,  ^  and,  F  =  ?  are  greater  than  unity,  with  Eq  the  peak  input  field.  Physically, 

this  requires  the  linear  absorption  to  be  large  on  a  wavelength  scale,  and  the  incident  field  to  be  strong 
enough  to  saturate  the  absorption.  In  the  case  of  ultrashort  pulses,  we  also  require  that  the  incident 
pulse  duration  tp  is  greater  than  the  polarization  decay  time  but  less  than  the  population  decay  time, 
ti  >  ip  >  t2  [1].  To  meet  these  conditions  for  the  sub- 100  fs  pulses  to  which  we  were  restricted  by 
computational  resources,  we  adopted  the  following  medium  parameters  for  illustrative  purposes:  ii  =  0.5 
ns,  t2  =  10  fs,  CO  =  U21  =  2  X  10^®  rads“^  (A  —  942nm),  p  —  4x  10“^®  Cm  and  No  =  2  x  10^°  cm”^.  For 
these  values  a  normalised  field  strength  of  F  =  1  corresponds  to  an  electric  field  strength  of  Eq  =  1.15  X  10® 
V/m.  The  linear  complex  refractive  index  of  the  absorbing  medium  is  given  by,  =  1  +  ?  where 

^  =  (^21  —  ^)/l2  is  the  normalised  detuning.  For  the  assumed  parameter  values  =  3.8,  and  has  a 
substantial  imaginary  contribution  so  that  the  linear  skin  effect  is  to  be  expected  at  the  absorber  boundary. 
The  nonlinear  Maxwell’s  e([uations  were  numerically  integrated  using  Yee’s  second-order  FD-TD  scheme 
for  the  field  updates  and  a  fourth-order  Runge-Kutta  method  for  the  medium  updates.  On  the  driving 
face  (z  =  0)  of  the  computational  domain,  the  incident  pulse  was  defined  to  be  Gaussian-shaped  in  space 
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and  time.  Periodic  boundary  conditions  were  imposed  at  the  transverse  boundaries  (ar  =  ±1^/2)  which 
were  placed  far  enough  away  for  the  optical  intensity  to  be  reduced  by  more  than  50dB  compared  with 
pulse  center.  The  medium  was  initialised  uniformly  with  p^i  =  0  and  n  =  1  for  z 

Figures  1  and  2  show  results  from  a  pair  of  sample  calculations  for  incident  pulses  with  peak  fields  of  (a) 
Eo  =  7  X  10^  V/m,  and  (b)  JSq  =  7  X  10®  V/m.  The  computational  domain  was  132  x  113  pm  {L^  X 
discretised  onto  a  3000  X  250  numerical  grid,  with  the  saturable  absorber  interface  situated  at  Za  =  125 
pm.  The  incident  pulse  was  80fs  long  (FWHM  of  the  electric  field),  with  a  beam  waist  of  Wq  =  20  pm 
(1/e  half- width  of  the  field),  and  was  initialised  such  that  its  peak  entered  the  computational  domain  at 
z  =  0  when  t  =  120/s.  The  temporal  snapshots  are  taken  after  reflection  from  the  interface  (at  t  =  770 
ps)  when  the  pulses  are  travelling  back  towards  z  =  0.  In  Fig.  1  only  the  extracted  field  envelopes  are 
shown,  since  the  carrier  is  not  easily  discernible.  The  corresponding  spectra  in  Fig.  2  were  obtained  from 
the  full  electric  field  data. 

Figures  1(a)  and  2(a)  show  that  for  low  incident  pulse  intensities  the  self-reflected  pulse,  though  diminished 
in  amplitude,  is  unchanged  spatially  or  spectrally.  The  pulse  energy  is  reduced  during  the  absorptive 
reflection,  and  the  calculated  reflectivity  of  24%  compares  favourably  with  the  predicted  steady-state  linear 
value  of  22%  from  Fresnel’s  law  and  the  expression  for  the  complex  refractive  index  above.  The  pulse  waist 
remains  unchanged,  evolving  almost  imperceptibly  on  the  time  scale  of  the  computation,  because  of  its 
long  diffraction  length  {Lo  =  ttwI/Xq  >  1000  pm).  In  this  case,  therefore,  since  the  pulse  is  too  weak 
to  significantly  saturate  the  absorber,  the  physical  properties  of  the  interface  remain  unmodified  and  the 
pulse  is  partially  reflected,  according  to  the  linear  skin  effect,  as  if  from  a  flat  mirror. 

In  contrast.  Figs.  1(b)  and  2(b)  show  that  the  high  intensity  pulse  undergoes  significant  nonlinear  spatial 
and  spectral  reshaping  at  the  interface.  Figure  1(b)  shows  the  pulse  profile  after  it  has  propagated  away 
from  the  interface,  close  to  its  focal  point,  where  the  pulse  waist  (at  maximum  Ey)  is  reduced  to  2.2  /xm, 
a  considerable  reduction  in  spot  size  in  comparison  to  the  linear  case  in  Fig.  1(a).  At  this  point,  the  peak 
electric  field  is  actually  fractionally  greater  than  the  peak  incident  field,  despite  the  decreased  reflectivity 
of  18%.  The  reflectivity  is  reduced  compared  to  the  linear  case  because  the  high  power  pulse  reduces  the 
medium  absorption  and  as  a  result  the  pulse  penetrates  further  into  the  interface. 

The  spectral  broadening  that  results  from  the  nonlinear  saturation  of  the  absorber  is  clearly  evident  in 
Fig.  2(b)  in  comparison  to  Fig.  2(a).  The  longitudinal  spectral  broadening  (in  Kz)  and  red-shifting  of  the 
reflected  field  are  due  to  the  self-reflection  of  the  incident  pulse  from  an  absorption  front  which  initially 
accelarates  and  then  deccelarates  into  the  saturable  medium,  which  in  turn  causes  the  reflected  light  to  be 
red-shifted  and  chirped. 

The  formation  of  a  curved  focusing  reflection  front  is  shown  in  Fig.  3,  where  the  spatial  distribution 
of  the  population  difference,  n,  is  plotted  after  its  creation  by  the  leading  edge  of  the  pulse.  Once  this 
curved  mirror  is  formed,  the  rest  of  the  pulse  is  reflected  towards  the  x  =  0  axis,  leading  to  the  tight  focus 
seen  in  Fig.  1(a).  The  spectral  width  of  the  broadest  feature  in  Fig.  2(b)  accurately  reflects  the  inverse 
incident  pulsewidth.  Furthermore,  we  note  the  signature  of  the  ordinary  (unfocused)  plane- wave  reflection 
at  the  leading  edge  of  the  ptdse  in  Fig.  1(a),  on  a  timescale  of  approximately  ^2?  during  which  the  mirror 
formation  takes  place. 
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Figure  Captions 

Fig.  1  Surface  plots  of  the  calculated  field  envelope  profiles  at  t=770fs  for  peak  input  fields  of  (a)  Eq  = 
7  X  10®  V/m  and  (b)  Eq  =  7  x  10*  V/m,  after  reflection  from  the  interface. 

Fig.  2  Contour  plots  of  the  spatial  Fourier  transforms  of  the  field  distributions  in  Fig.  1. 

Fig.  3  The  population  difference  of  the  saturable  absorber  following  self-reflection  of  the  Eq  =  7  X  10®  V/m 
(high  power)  pulse  in  Fig.  1(b)  at  t  =  770  fs.  The  contour  at  n  =  0.5  is  superimposed. 
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Figure  1 
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ID  and  2D  Femtosecond  Light  Bullets  and  Bubbles 
S,  Glasgow,  J.  Xin,  N.  Ercolani  and  J.V.  Moloney 

The  existence  of  localized  (non-coUapsing)  femtosecond  duration  optical  pulses  with  very  high  local  inten¬ 
sities  but  very  low  energies,  presents  a  fascinating  opportunity  to  revisit  the  optical  breakdown  problem 
which  has  plagued  high  power  laser  design  for  decades.  In  many  instances,  the  envelope  approximations 
to  the  vector  Maxwell  equation  predict  coUapse  to  a  singularity  in  finite  time.  An  open  issue  is  whether 
higher-order  corrections  to  envelope  approximations  can  introduce  regularization  of  blow-up  or  whether 
one  has  to  resort  to  a  fuU-fledged  Maxwell  description.  We  have  been  studying  these  phenomena  in  ID 
and  2D  with  the  view  to  establishing  a  connection  between  these  distinct  levels  of  treatment.  In  ID  we 
have  explored  localization  with  an  instantaneous  cubic  nonlinear  term  and  the  analog  of  critical  collapse 
with  an  instantaneous  qtiintic  nonlinear  term. 

Maxwell’s  equations  have  the  following  form: 


X  V  X  ^ 

(1) 

+  ^o'P  =  (^5  “ 

(2) 

11 

1 

(3) 

which  allows  for  both  a  linear  and  nonlinear  polarizations  within  a  simple  quantum  oscillator  model.  This 
model  was  initially  used  to  revisit  many  well-established  nonlinear  optical  interactions  but  now  occurring  on 
femtosecond  timescales.  The  validity  of  traditional  envelope  approximations  becomes  an  issue  as  the  pulse 
gets  progressively  shorter.  Moreover,  as  the  pulse  shortens,  its  energy  decreases  and  mucher  higher  peak 
intensities  can  be  achieved  before  material  breakdown  becomes  an  issue.  The  need  for  a  full  reappraisal  of 
light-material  coupling  on  such  ultra-short  timescales  becomes  essential. 

Our  recent  work  showed  that  the  assumption  of  an  instantaneous  nonlinear  optical  response  would  lead 
to  shocking  of  femtosecond  pulses.  We  have  now  adopted  the  more  realistic  quantum  oscillator  model 
above,  which  includes  linear  and  nonlinear  dispersion,  and  from  which  all  of  the  known  classical  envelope 
descriptions  can  be  obtained  via  singular  perturbation  expansions.  We  have  confirmed  the  existence  of  a 
2D  TE  light  bullet  and  have  evidence  to  show  that  a  recently  re-discovered  ID  unipolar  (half-cycle)  light 
bubble  is  unstable  in  higher  space  dimensions.  The  ID  version  of  this  latter  object  has  been  shown  to  be 
integrable  by  us.  Figure  1  shows  numerically  computed  2D  TE  light  bullet  (a)  and  bubble  (b).  The  latter 
is  decaying  slowly. 


E(x,2.  t  =  30  fs) 
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Large  Scale  Computational  Nonlinear  Optics 

Maxwell  semiconductor  Bloch  equations 
R.  Indik  and  J.V.  Moloney 

To  model  the  behavior  of  nltrashort  pulses  in  semiconductor  media,  we  use  the  Maxwell- Semiconductor 
Bloch  equations.  As  mentioned  earlier  in  this  report,  these  equations  describe  how  the  electric  field  E 
evolves  in  the  traveling  frame  of  reference  (  =  z  -  tcfn.  In  practice,  we  assume  isotropy  in  q,  and 
motivated  by  device  geometries  assume  strong  guiding  in  one  of  the  transverse  spatial  dimentions.  This 
reduces  the  problem  to  calculation  of  the  dynamics  of  integro- differential  equation  in  3  dimensions. 


/p  2. 

2ikodx^  koV  4  ^ 
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This  evolution  is  given  in  terms  of  the  momentum  resolved  polarization  Pg.  The  evolution  of  this  polar¬ 
ization  and  the  electron/hole  distribution  functions  n®/'*  in  time  t  is  given  by 
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with  the  renormalized  Rabi  frequency 

ng  =  -^  +  —  I  q'Vg,g.Pg,q^dq 

and  renormalized  energy  dispersion 

^  ~  W  Ig> 

Vg^gi  is  the  angle  averaged  coulumb  interaction.  The  details  of  the  collision  terms  are  described  in  an  earlier 
section  of  this  report.  They  have  little  impact  on  the  methodology  of  the  numerical  implementation. 
These  are  differential  equations  in  t,  ^  and  x,  and  integral  in  q,  the  magnitude  of  the  momentum.  The 
integrals  in  this  calculation  lead  to  a  tight  coupling  of  the  equation  across  the  momentum  variable  q.  This 
is  precisely  as  a  result  of  the  inclusion  of  the  Coulumb  interaction,  which  enters  in  these  equations  through 
the  band-gap  renormalization  Ag  and  the  generalized  Rabi  frequency  OmegUg. 

Numerical  integration  of  the  MaxweU-Semiconductor-Bloch  (MSB)  equations  quickly  becomes  a  large  prob¬ 
lem.  Even  in  the  plane-wave  case,  where  a  strong  transverse  guiding  assumption  in  x  reduces  the  system 
from  dynamics  in  three  dimensions  to  two  dimensions,  the  parabolic  band  structure  leads  to  the  necessity 
of  integrating  the  behavior  of  tightly  coupled  oscillators  whose  frequencies  vary  over  a  large  fraction  of  the 
underlying  optical  frequency,  each  of  which  has  a  decay  rate  of  the  order  of  50  fs.  In  practice,  this  forces 
explicit  step  sizes  to  be  of  order  0.5  fs.  If  we  are  primarily  interested  in  the  behavior  of  very  thin  slabs 
of  semiconductor,  the  comi)utations  are  not  too  onerous.  However  propagation  through  thick  samples  is 
another  matter,  and  when  the  transverse  guiding  assumption  is  relaxed,  and  the  3  dimensional  problem 
with  propagation  is  considered,  the  calculation  becomes  quite  challenging. 

We  have  implemented  a  code  that  can  follow  the  three  dimensional  dynamics  of  a  pulse  in  a  slab  amplifier 
without  transverse  confinement.  This  code  has  been  developed  using  PVM  and  run  on  the  IBM  SP2  at  the 
Mauii  High  Performance  Computing  Center  (MHPCC),  on  the  Thinking  Machines  CM5  at  the  Army  High 
Performance  Computing  Center  at  the  Minnesota  Supercomputing  Center,  as  well  as  on  a  network  of  Sun 
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Spare-stations,  Silicon  Graphics  Indy  work  stations  and  a  four  processor  Silicon  Graphics  Power  Challenge 
at  the  University  of  Arizona.  The  code  takes  advantage  of  the  natural  parallelism  of  the  problem  in  the 
transverse  dimension  to  achieve  very  efficient  parallel  execution.  The  same  code  has  been  used  to  make 
exploration  in  parameter  space  for  the  plane  wave  problem.  The  problem  is  discretized  in  momentum, 
time  and  space,  and  is  solved  using  the  method  of  lines.  The  time  equations  are  integrated  using  fourth 
order  Runge-Kutta,  while  t].ie  spatial  propagation  is  done  using  second  order  central  differences  and  second 
order  Runge-Kutta.  A  typical  long  propagation  computation  including  transverse  structure  uses  2.5-5  epu 
hours  on  64  nodes  of  the  IBM  SP2  at  MHPCC. 

The  structure  of  the  code  is  illustrated  in  figure  1.  We  use  a  domain  decomposition,  to  parcel  the  work 
of  calculation  over  many  processors.  Each  processor  is  responsible  for  calculating  the  propagation  of  the 
E{x,t,^)  at  a  point  x  =  x,-.  The  only  place  in  the  calculation  which  requires  that  processors  exchange 
information  (other  than  the  initialization  and  for  output)  is  in  the  calculation  of  the  diffraction  term 
d^E/dx^.  Since  we  are  using  second  order  central  differences,  this  requires  nearest  neighbor  exchange  of 
the  values  of  E.  This  is  a  very  small  amount  of  information  that  must  be  communicated  compared  to  the 
amount  of  data  in  the  q  resolved  polarization  P,  and  densities  Not  only  is  little  data  shared  between 
processors,  that  data  is  only  shared  relatively  infrequently.  The  vast  bulk  of  the  time  that  is  required  for 
the  computation  is  spent  in  integrating  the  local  material  equation  in  time  t.  The  values  of  the  field  E 
can  be  sent  to  nearest  neighbors  before  the  calculation  of  the  material  response  is  started,  so  that 

the  propagation  of  P  in  ^  (requiring  the  values  of  the  neighboring  E  values)  is  done  after  the  material 
response  has  been  computed.  Thus  even  though  communication  between  processors  can  be  quite  slow, 
the  computation  need  not  wait.  In  practice,  more  than  one  tarnsverse  value  may  be  executing  on  one 
processor,  as  we  take  advantage  of  the  virtual  machines  that  PVM  provides.  This  allows  us  to  run  the 
same  code  on  our  local  network  as  we  run  on  the  remote  sites. 

Our  experience  using  the  The  CM-5  at  the  AHPCC  of  the  MSC  (Minnesota  Supercomputing  Center)  and 
the  SP2  at  the  MHPCC  (Mauii  High  performace  Computing  Center)  has  led  us  to  conclude  that  for  the  very 
coarse  grained  parallelism  in  our  code,  the  SP2  in  dedicated  mode  is  nearly  ideal.  The  individual  machines 
(processors)  execute  the  code  independently  and  quite  efficiently,  and  the  communication  overhead  is 
minimal.  As  the  usage  at  MHPCC  has  increased,  the  performance  of  our  code  has  degraded  dramatically. 
This  is  due  to  the  uneven  load  across  the  many  machines.  At  any  given  time,  our  code  executes  at  the 
speed  of  the  slowest  of  the  machines.  In  addition,  although  dedicated  machine  time  has  been  available,  the 
waiting  time  to  get  that  time  has  been  quite  long.  We  are  hopeful  that  the  implementation  LoadLeveler  at 
MHPCC  wiU  improve  this  situation  dramatically.  This  should  be  installed  in  mid  March  1995.  We  have 
found  that  the  user  support  services  at  both  MSC  and  particularly  at  MHPCC  have  been  excellent. 
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Slaves  do  all  of  the  work,  and  communicate  among  themselves  to  calculate  diffraction 
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In  the  propagation  of  the  light  field  E(t),  diffraction  couples 
nearest  neighbors. 
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Large  Scale  Computational  Nonlinear  Optics  -  Maxwell’s  Equations 
R.G.  Flesch,  and  J.V.  Moloney 

Solution  of  Maxwell’s  equations  in  three  spatial  dimensions  plus  time  places  such  memory,  speed  and  disk 
storage  demands  on  the  systems  as  to  put  it  is  the  class  of  Grand  Challenge  projects  like  lattice  gauge 
theory.  Resolutions  of  only  20  points  per  wavelength  in  the  longitudinal  direction  can  lead  to  executables 
as  large  a.s  10  Gbytes.  Clearly  the  only  machines  which  have  such  resources  are  the  largest  supercomputers 
such  as  the  CM5  at  the  AHPRC  and  the  SP2  in  Maui.  That  is  not  to  say  that  one  cannot  make  good 
progress  with  more  modest  resources.  Our  strategy  has  been  to  develop  the  codes  on  our  Power  Challenge 
and  do  preliminary  parameter  searches  with  smaller  grid  resolution.  Once  an  interesting  regime  has  been 
uncovered  we  can  test  and  further  probe  the  phenomena  on  the  larger  machines  at  increased  resolution  to 
ensure  accuracy.  A  sketch  of  the  computational  domain  illustrating  one  means  of  visualization  in  movie 
format  is  shown  in  the  figure. 

Algorithmically,  the  solution  process  cleanly  breaks  into  two  steps,  first  the  Maxwell  field  updates  followed 
by  computation  of  the  material  response.  The  first  of  these  two  steps  requires  a  great  deal  of  commu¬ 
nication  while  the  second  in  completely  local.  Both  steps  require  roughly  the  same  amount  of  time,  the 
communication  in  the  Maxwell  steps  is  balanced  by  the  more  intense  number  of  actual  flops  in  the  material. 
As  is  the  case  in  any  multi-processor  machine,  each  processor  is  allocated  an  equal  fraction  of  grid  points. 
To  minimize  interprocessor  communication,  it  is  essential  to  minimize  the  surface  area  to  volume  ratio, 
something  which  is  the  default  data  layout  on  the  CM5.  The  high  interprocessor  bandwidth  and  relatively 
low  latency  make  the  CM5  the  optimal  choice  we  have  examined  to  date.  We  have  yet  to  port  the  code 
to  the  SP2  as  the  PVM  paradigm  will  require  complete  recoding.  The  success  of  the  port  will  depend  on 
how  well  the  high  speed  switch  on  the  SP2  handles  the  large  amount  of  communication.  As  the  second 
step  requires  no  communication,  the  CM5,  C90  and  Power  Challenge  all  parallelize  that  section  of  the  code 
very  well. 

One  very  important  numerical  aspect  of  the  problem  which  we  had  to  address  is  the  presence  of  numerical 
dispersion  in  the  finite  differencing  schemes  used  in  the  Maxwell  step.  There  have  recently  been  a  significant 
number  of  publications  on  this  issue.  We  have  run  our  own  tests  using  a  model  system  in  which  there  is 
no  physical  dispersion.  We  have  observed  the  phase  shifts  reported  in  the  literature  using  both  the  (2,2) 
and  (2,4)  numerical  schemes  and  have  been  able  to  control  it  by  using  a  larger  number  of  grid  points. 
However  we  have  found  that  the  actual  shifts  in  the  position  of  the  packets  we  propagate  to  be  very 
small.  For  example  have  found  no  real  difference  in  the  carrier  shock  times  if  we  use  200  or  1000  points 
per  wavelength.  This  is  in  great  part  due  to  the  fact  that  we  use  ultrashort  pulses  and  propagate  them 
hundreds  or  a  thousand  wavelengths.  The  process  is  cumulative  and  would  introduce  more  errors  for  longer 
propagations.  Equally  important  is  that  for  realistic  material  models  the  physical  dispersion  far  exceeds 
the  numerical. 

Another  aspect  which  is  essential  to  a  successful  implementation  of  this  research  program  is  the  visualization 
of  the  large  date  sets  generated.  Visualization  of  a  three  dimensional  three  component  vector  field  has 
always  been  a  difficult  task  and  one  must  normally  resort  viewing  cross  sections.  One  software  tool  which 
greatly  eases  this  task  is  AVS,  enabling  one  to  load  in  an  entire  3D  data  set  and  view  any  desired  slice. 
But  often  this  also  does  not  suffice  because  due  to  the  rapid  longitudinal  variation  due  to  the  carrier  can 
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obstruct  field  values.  To  avoid  this  we  have  written  some  GL  code  for  our  Silicon  Graphics  machines.  With 
this  code  we  are  able  to  rotate  and  animate  two  dimensional  surfaces  of  the  field  avoiding  the  obstruction 
problems.  The  hardware  graphics  present  in  the  Indigo  2  allow  us  to  do  these  rotations  in  real  time. 
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Relaxation  processes  in  semiconductor  light  amplifiers:  the  role  of  charge  carrier  cor¬ 
relations 

Y.  Lvov,  R.  Binder  and  J.V.  Moloney 

One  of  the  most  important  aspects  of  a  microscopic  approach  to  semiconductor  lasers  and  semiconductor 
light  amplifiers  is  the  problem  of  charge  carrier  relaxation  and  thermalization  processes.  Being  far  away 
from  thermal  equilibrium,  the  coupled  semiconductor-light  system  is  generally  not  well  described  in  terms 
of  equilibrium  concepts.  Many  theories  use  the  concepts  of  quasi-thermal  equilibrium  amended  by  phe¬ 
nomenological  thermalization  terms.  In  the  quasi-thermal  equilibrium  approach,  the  conduction  band  elec¬ 
trons  and  valence  band  holes  are  assumed  to  be  independently  in  equilibrium,  and  the  thermalization  terms 
account  for  carrier- carrier  and  carrier-phonon  scattering  which  tends  to  establish  the  quasi-equlibrium  if, 
for  example  due  to  the  presence  to  the  light  field,  the  system  is  not  in  the  quasi-equilibrium  state.  Since  the 
incoherent  scattering  processes  crucially  affect  the  optical  gain  dynamics  of  the  systems  (in  other  words, 
determine  the  filling  of  energy  states  which  contributes  to  stimulated  emission),  a  microscopic  treatment 
of  these  terms  is  clearly  desirable. 

In  recent  years,  several  research  groups  have  already  implemented  basic  versions  of  microscopic  scattering 
models  in  their  analysis  of  optical  properties  of  semconductors.  Most  of  these  models  are  based  on  the 
combination  of  three  approximations:  the  screened-Hatree-Fock  approximation  (i.e.,  the  scattering  prob¬ 
ability  being  proportional  to  the  square  of  the  screened  Coulomb  interaction  matrix  element,  and  being 
solely  determined  by  carrier  distribution  functions,  not  carrier  correlation  functions),  the  Markov  approx¬ 
imation,  and  the  completed  scattering  approximation  (i.e.,  the  scattering  contains  an  energy- conserving 
delta-function).  These  are  typical  features  of  the  so-caUed  Kadanoff-Baym  theory,  which  was  developed  to 
describe  systems  very  close  to  equilibrium. 

In  systems  which  deviate  considerably  from  equilibrium,  such  as  lasers  operating  far  above  threshold,  lasers 
with  fast  external  modulations  and  light-amplifiers  amplifying  sub-ps  pulses,  these  approximations  become 
questionable.  Also,  novel  laser  systems  like  ZnSe-based  lasers  cannot  simply  be  described  on  the  basis  of 
commonly  used  weakly-correlated  plasma  theories.  Although  the  lasing  processes  in  these  structures  are 
stiU  being  investigated  (with  special  emphasis  on  excitonic,  biexcitonic,  and  carrier-phonon  processes),  it 
seems  appropriate  at  the  present  time  to  amend  existing  theories  and  include  charge-carrier  correlations 
beyond  the  screened-Hartree-Fock  approximation. 

In  the  present  project  we  are  deriving  equations  of  motion  for  carrier-correlation  functions,  which  ultimately 
should  allow  us  to  study  relaxation  mechanisms  in  systems  with  strong  Coulomb  correlations. 

Due  to  the  nature  of  the  quantum  mechanical  Coulomb  interaction  Hamiltonian,  the  time  evolution  of  the 
distribution  functions  (two  operator  expectation  values)  n*  =<  afak  >  are  determined  by  four  operator 
expectation  values  ^4(^1,  A;2,  ^3,  A;4)  =<  >.  The  time  evolution  of  the  four  operator  expecta¬ 

tion  values  is  defined  by  six  operator  expectation  values,  and  so  on.  This  infinite  hierarchy  of  operator 
equations  is  known  as  BBGKY  hierarchy.  In  order  to  obtain  finite  sets  of  equations  one  has  to  decouple  the 
hierarchy  with,  for  example,  the  help  of  Green’s  function  formalisms.  This  is  done  in  the  Kadanoff-Baym 
approach.  Using  this  formalism  one  can  obtain  Boltzmann  kinetic  equations,  which  within  the  certain 
limits  describe  the  time  evolution  of  two  operator  expectation  values.  The  Green’s  function  approach,  al¬ 
though  being  a  very  practical  approach,  has  two  disadvantages:  the  two-time  structure  of  the  one-particle 
Green’s  functions  and  the  lack  of  precise  and  systematic  error  studies  of  given  approximations  (i.e.,  without 
going  beyond  the  screened-Hartree  Fock  approximation  no  definite  statements  about  its  accuracy  can  be 
made). 

We  propose  an  alternative  approach  for  studying  carrier-carrier  collisions  and  correlations  which  combines 
the  method  of  nonlinear  optical  suseptibilities  for  Coulomb  systems,  recently  developed  by  Axt  and  Stahl) 
with  the  method  of  multiple  scales.  The  combination  of  these  approaches  should  enable  us  to  obtain 
generalized  Boltzmann  equations  including  charge- carrier  correlation  functions.  As  this  project  is  based 
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on  a  highly  demanding  and  complex  theory,  we  are  planning  the  following  a  two  step  approach. 

In  the  first  step,  we  neglect  aU  contributions  of  the  optical  field  and  study  correlations  functions  in  a 
certain  order  of  the  interaction  potential  (which,  at  this  point,  we  assume  to  be  small).  We  can  then 
apply  the  multiple  scales  method  to  obtain  an  equation  of  motion  for  the  four  operator  expectation  value 
SN4  =  N4  -  n  X  n,  where  n  x  n  denotes  the  corresponding  Hartree-fock  factorized  form  of  N4.  Initial 
calculations  indicate  a  strong  formal  similarity  of  the  generalized  Boltzmann  equation  with  the  conventional 
Boltzmann  equation. 

We  are  also  interested  in  studying  stationary  solutions  of  the  systems  described  by  the  above  mentioned 
relaxation  equations  and  additional  (more  or  less  phenomenological)  carrier  source  and  sink  terms.  These 
questions  are  related  to  the  well-known  Kolmogorov  solutions  in  conventional,  classical  scattering  investi¬ 
gations. 

After  acquiring  a  certain  familiarity  with  the  temporal  evolution  of  correlation  functions  in  weakly  inter¬ 
acting  systems  (i.e.  sytems  which  allow  for  a  classification  of  solutions  corresponding  to  the  order  of  the 
interaction  potential),  we  plan  to  generalize  the  approach  and  include  the  optical  polarization,  similar  the 
the  Axt-Stahl  approach.  Of  course,  in  contrast  to  their  investigation,  we  have  to  go  to  higher  than  third 
order  in  the  optical  field  to  properly  describe  relaxation  effects. 


46 


Feedback  and  Multimode  Behavior  of  a  Fabry-Perot  Laser:  A  Travelling  Wave  Model 
M.  Homar,  J.V.  Moloney  and  M.  San  Miguel 

Summary 

We  report  the  results  of  a  Lumerical  study  of  multimode  behavior  of  a  Fabry-Perot  laser[l].  The  model  is 
based  on  travelling-wave  equations  for  the  slowly  varying  amplitudes  of  the  counterpropagating  waves  in 
the  cavity,  coupled  to  equations  for  spatially  dependent  population  inversion  and  polarization  of  a  two-level 
active  medium.  Variations  in  polarization  and  population  inversion  over  wavelength  scales  are  taken  into 
account  by  means  of  an  expansion  in  a  Fourier  series.  Results  are  given  for  typical  semiconductor  Iciser 
parameters.  Spatially  distributed  spontaneous  emission  noise  and  carrier  diffusion  are  taken  into  account. 
The  competing  role  of  Spectral  Hole  Burning  (SHB),  spontaneous  emission  noise  and  carrier  diffusion 
in  determining  multimode  behavior  is  elucidated;  With  no  carrier  diffusion,  spontaneous  emission  noise 
excites  a  large  number  of  modes  close  to  threshold,  while  SHB  leads  to  a  fixed  number  of  significant 
lasing  modes  well  above  threshold.  Carrier  diffusion  restores  dominant  single-mode  emission  well  above 
threshold.  We  have  also  studied  the  effects  of  optical  feedback  and  opportunities  for  mode  selection  with 
short  external  cavities:  For  a  cavity  of  Lgxt  =  75;um  a  chosen  single  mode  can  be  selected  for  an  external 
field  mirror  reflectivity  of  1%.  For  a  reflectivity  of  40%  two  groups  of  intracavity  modes  separated  by  the 
external  cavity  mode  interspacing  are  selected.  For  a  cavity  of  L^xt  =  450^m  the  laser  can  be  forced  to 
lase  in  two  intracavity  modes  for  a  reflectivity  of  2%.  Mode  selection  is  not  found,  even  for  weak  feedback, 
for  external  cavities  longer  than  0.1cm. 


Traveling  wave  model 

The  interaction  of  the  electromagnetic  wave  with  the  medium  is  governed  by  the  Bloch  equations  for 
two-level  atoms. 
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In  the  equations  above  we  have  also  included  Langevin  noise  terms  ^(p){zyt)  as  polarization  sources.  Such 
spatially  distributed  noise  terms  model  independent  spontaneous  emission  process  in  different  points  of 
the  cavity.  They  are  taken  to  be  Gaussian  white  noise  in  space  and  time  with  zero  mean  and  correlations 
(^o(^5  0^o(^^^0)  “  “  f)6{z  “  z^).  Longitudinal  carrier  diffusion  is  also  taken  into  account. 

We  denote  by  ri  and  the  reflectivities  of  the  cavity  mirrors  placed  at  ^  =  0  and  z  =  respectively.  In 
addition  we  consider  an  external  mirror  of  reflectivity  Vs  placed  on  the  right  of  the  laser  cavity.  Thus,  our 
boundary  conditions  are, 

E'^(z  =  0^t)  =  riE’^{z  =  0^t) 

E~{z  =  L,t)  =  r2E^{z  =  L, /)  -  rlE^{z  =  L^t  -  r) 
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where  r  is  the  round  trip  time  of  the  field  in  the  external  cavity  of  length  LeXt.  Note  that  these  boundary 
conditions  take  into  account  multiple  reflections  on  the  external  mirror. 
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Circularly  symmetric  distributed  feedback  semiconductor  lasers 
Ewan  M.  Wright  and  K.  Kasunic 

We  are  currently  investigating  the  modal  properties  of  circularly  symmetric  distributed  feedback  (CSDFB) 
lasers  using  both  coupled-mode  and  beam  propagation  techniques.  CSDFB  lasers  offer  a  number  of  ad¬ 
vantages  over  conventional  surface  emitting  lasers  including  circularly  symmetric  output,  low- divergence 
output  beam,  and  the  ability  to  phase-lock  a  two-dimensional  array. 

Previous  work  on  CSDFB  lasers  concentrated  on  the  mode  properties  and  the  calculation  of  the  laser 
threshold  using  coupled-mode  theory.  We  have  extended  this  work  to  deal  explicitly  with  issues  related  to 
CSDFB  semiconductor  lasers  including  index  anti-guiding  associated  with  the  gain  profile,  azimuthal  mode 
selection,  and  the  laser  linewidth.  The  key  issue  we  have  addressed  is  that  of  the  azimuthal  mode  selection 
since  this  determines  whether  or  not  the  output  is  circularly  symmetric.  However,  as  one  increases  the 
azimuthal  mode  selection  we  find  that  the  laser  linewidth  also  increases,  meaning  that  as  a  result  more 
azimuthal  modes  may  actually  oscillate  close  to  threshold!  Our  initial  result  indicate  that  CSDFB  lasers 
may  typically  oscillate  on  several  azimuthal  modes,  and  hence  have  non-circular  output  beams,  close 
to  threshold.  Far  above  threshold  the  laser  linewidth  narrows  and  the  output  beam  may  then  become 
circularly  symmetric. 

To  extend  our  understanding  of  CSDFB  semiconductor  lasers  we  are  currently  developing  beam  propagation 
codes  which  also  incorporate  the  nonlinear  optical  response  of  the  semiconductor  medium.  This  wiU  allow 
us  to  address  the  problem  of  nonlinear  transverse  mode  selection  and  also  the  laser  stability  which  is  beyond 
the  previous  coupled  mode  theories. 

Our  most  recent  work  has  concentrated  on  the  above-threshold,  nonlinear  properties  of  CSDFB  lasers. 
Specific  issues  addressed  in  our  research  are  gain  saturation  [1],  index  saturation  [2],  and  solutions  to  the 
time-dependent  coupled-mode  equations  on  a  finite  interval. 

By  including  radially-vaxying  gain  saturation  and  uniformly-distributed  surface  emission  in  the  above¬ 
threshold  solution  of  the  coupled-mode  equations,  we  have  obtained  2  significant  results.  First,  there  exists 
an  optimum  coupling  strength  k  for  efficient  conversion  of  pump  power  into  useful  output  power.  From 
the  perspective  of  device  design,  these  results  allow  appropriate  selection  of  k  for  high-power  operation, 
alleviating  potential  problems  with  overheating.  Second,  for  relatively  large  k,  high-power  operation  of  the 
CSDFB  laser  wiU  result  in  multi-mode  operation,  a  consequence  of  an  above-threshold  mode  competition. 
Fortunately,  this  multi-mode  behavior  can  be  avoided  by  designing  these  lasers  with  the  optimum  k  for 
power  conversion. 

We  have  also  modeled  the  effects  of  gain/index  coupling,  by  including  index  saturation  via  the  Unewidth 
enhancement  factor  a.  The  results  here  show  that  even  for  large  k,  single-mode  behavior  is  possible  over 
a  Umited  (above-threshold)  power  range.  We  also  believe  our  results  help  explain  the  difficulties  found  in 
experimentaUy  obtaining  above-  threshold,  single-mode  output. 

FinaUy,  we  have  started  the  numerical  modeUng  of  the  time-dependent  coupled-mode  equations  describing 
self-pulsing  and  soliton  formation  in  nonlinear  periodic  structures.  To  date,  we  have  successfuUy  modeled 
self-  pulsing  in  conventional  DFB  structures.  Our  future  work  will  extend  these  results  to  the  CSDFB 
laser. 
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Abstract 

Quasi-equilibrium  excitation  dependent  optical  probe  spectra  of  II- VI  semi¬ 
conductor  quantum  wells  at  room  temperature  are  investigated  within  the 
framework  of  multi-band  semiconductor  Bloch  equations.  The  calculations 
include  correlation  effects  beyond  the  Hartree-Fock  level  which  describe  de¬ 
phasing,  interband  Coulomb  correlations  and  band-gap  renormalization  in 
second  Born  approximation.  In  addition  to  the  carrier-Coulomb  interaction, 
the  influence  of  carrier-phonon  scattering  and  inhomogeneous  broadening  is 
considered.  The  explicit  calculation  of  single  particle  properties  like  band 
structure  and  dipole  matrix  elements  using  k  •  p  theory  makes  it  possible  to 
investigate  various  II- VT  material  combinations.  Numerical  results  are  pre¬ 
sented  for  CdZnSe/ZnSe  and  CdZnSe/MgZnSSe  semiconductor  quantum-well 
systems.  _ 

Keywords:  A.  semiconductors  B.  optical  properties  C.  carrier-carrier  interaction 
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I.  INTRODUCTION 


Wide  band  gap  semiconductor  materials,  in  particular  II- VI  compounds,  are  promising 
candidates  for  optical  device  application  in  the  blue-green  wavelength  region.  Especially 
for  device  optimization,  a  detailed  understanding  of  the  electronic  interaction  processes 
and  their  influence  on  optical  gain  spectra  is  desirable.  Several  recently  investigated  II- VI 
heterostructures  are  composed  of  quantum  layers  with  Znj;Cdi_a;Se  active  material  within 
ZnSe  barriers  [1-4].  Hence,  the  theoretical  description  of  the  optical  properties  of  II- VI 
heterostructures  requires  not  only  the  inclusion  of  the  relevant  many-body  processes  in  a 
highly  excited  semiconductor  but  also  the  particular  composition  and  geometry  dependen^  -t 
band  structure  properties. 

In  this  paper  we  theoretically  investigate  the  optical  probe  spectra  of  II- VI  heterostruc¬ 
tures.  Concentrating  on  room  temperature  properties  of  structures  with  weak  interface 
roughness,  the  optical  gain  can  be  assumed  to  result  from  electrons  which,  after  their  inco¬ 
herent  injection,  have  relaxed  into  the  available  electronic  states  at  the  bottom  of  the  band 
whereas  higher  bound  states,  such  as  biexcitons  [5],  are  of  minor  importance.  In  our  model, 

v/* 

the  interaction  of  the  light  field  with  the  inverted  material  is  treated  semiclassically  by  cal¬ 
culating  the  absorption  in  the  framework  of  the  multi-band  semiconductor  Bloch  equations. 
The  observable  in  a  typical  experiment,  where  the  probe  light  propagates  in  the  plane  of  the 
quantum-well  heterostructure,  is  the  absorption  coefficient. 


where  rib/the  refractive  index,  E{u)  =  J  dt  E{t)  is  the  Fourier  transform  of  the  applied 
probe  pulse  E{t)  and  P{u))  is  the  corresponding  Fourier  transform  of  the  induced  polarization 
density  P{t).  Note  that  for  a  weak  probe  field  the  induced  polarization  P{lo)  is  directly 
proportional  to  E{u})  so  that  a(a;)  is  independent  of  the  probe-pulse  parameters. 

When  light  propagation  in  the  quantum-well  plane  is  experimentally  investigated,  the 
field  overlapjl^  with  the  quantum  well(s)  and  the  barrier  layers.  To  account  for  the  resulting 


2 


reduced  effective  absorption  or  gain,  we  have  introduced  in  Eq.  (1)  the  confinement  factor 
r  =  0.003/nm. 


II.  EQUATIONS  OF  MOTION 


In  this  section,  we  calculate  the  polarization  density  P{t)  from  a  microsocopic  theory 
including  band-structure  effects  as  well  as  many  body  interactions.  We  expand  the  polar¬ 
ization  density  in  a  quantum-well  Bloch  basis. 


^W  =  4EEA^r>r'(i)+c.c., 


(2) 


Pyl/'  k 


where  k  is  the  in-plane  carrier  momentum,  i/  =  A,  n  contains  the  band  index  X  =  e,h  and 
the  sub-band  index  n  and  A  is  the  active  area  of  the  quantum  well.  In  Eq.  (2),  we  consider 
all  optically  allowed  single  valence  band  to  conduction  band  (|z/)  — >  \v))  transitions  with 


the  interband  dipole  matrix  elements  .  The  polarizations  for  the  different  bands  obey 


the  multi-band  semiconductor  Bloch  equations  [6], 


- 1  [i - /k  -  /k ] + 


n 


dt 


(3) 


where  P^''  (t)  is  the  interband  polarization  between  the  bands  u  =  e,n  and  i/  =  h,n  and 

/k’”)  /k’”  carrier  distribution  functions  for  electrons  and  holes  in  the  sub-band  n, 

respectively.  For  a  weak  probe  beam,  the  polarization  P^''^  {t)  remains  linear  with  respect  to 

the  probe  field  E{t)  and  the  probe  field  induced  changes  of  the  carrier  occupation  can  be 

neglected.  Hence  the  probe  spectrum  reflects  the  excitation  of  the  system  which  is  described 

by  quasi-equilibrium  Fermi-Dirac  distribution  functions  with  fixed  carrier  temperature 

(T=300K)  and  a  chemical  potential  determined  by  the  total  density  in  the  sample. 

i_  __ — - - - - 

In  Eq.  (2),  the  many-body  effects  can  be  divided  into  Hartee-Fock  (mean  field)  and 

correlation  contributions.  The  mean  field  corrections  lead  to  renormalized  single-particle 

energies. 


=  e„(k)  -  €..(k)  -  Y,  Mk-ql  te  +  f^) ,  (4) 

q?tk 
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and  to  a  renormalized  Rabi  energy, 

flk  w v<.i  (5) 

q#k 

with  the  free-carrier  energies  e„(k)  and  the  Rabi  energy  of  the  probe  field  E.  The 
particular  properties  of  the  heterostructure,  such  as  quantum-well  thickness  and  material 
composition,  determine  the  dipole  matrix  transition  elements  as  well  as  the  bandstruc- 
ture  Cv(k).  In  our  approach,  these  quantities  are  calculated  from  a  diagonalization  of  the 
Luttinger  Hamiltonian  using  a  4x4  k*  p-theory  within  the  envelope  approximation  [6]. 

The  quantum-well  matrix  elements  of  the  bare  Coulomb  potential,  which  couple  various 
carrier  states,  have  the  general  form 

V^'^''‘'’‘0cuk2,q)  =  ^  fdzdz'  (6) 

X  e-'I'-'l  +  q,/)  -  q,z). 

Because  of  the  band  mixing,  the  confinement  functions  of  the  band  Ui  depend  on 

the  in-plane  carrier  momentum  ki.  This  band  mixing  is  additionally  complicated  by  the  fact 
that  the  influence  of  barrier  states  on  the  top  valence  band  states  cannot  be  neglected.  The 
dependence  of  the  quantum-well  Coulomb  matrix  elements  on  the  band  indices  i/,-  and  the 
corresponding  in-plane  carrier  momenta  ki  considerably  complicates  the  coupling  of  various 
bands  in  Eqs.  (3)-(5)  as  well  as  the  evaluation  of  screening.  It  would  be  desirable  to  use 
Eq.  (6)  without  further  approximations  within  the  many-body  problem.  Then  one  would 
take  into  account  that  the  Coulomb  interaction  explicitly  depends  on  the  momenta  of  the 
contributing  carriers  ki  and  not  only  on  the  transition  momentum  q.  In  this  paper,  we  study 
the  influence  of  band  mixing  under  the  assumption  that  the  various  envelope  functions 
^i/.(k,>2;)  in  Eq.  (6)  are  approximated  by  an  effective  (momentum  independent)  function 
The  resulting  confinement  Coulomb  potential 

H  =  —  / dzdz'  Mz)f  e-’l-*'!  |eo{/)r  (7) 

will  be  used  in  Eqs.  (4),  (5)  and  below. 
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The  mean-field  contributions  contain  the  lowest-order  Coulomb  effects  which  lead  to 
excitonic  resonances  and  a  low-density  band-gap  renormalization.  However,  under  high 
excitation  conditions,  correlation  contributions  [7-11]  result  in  strong  modifications  of  the 
Hartree-Fock  terms.  The  correlation  terms  can  be  divided  into  diagonal  (d)  and  non-diagonal 
(nd)  contributions: 

=  -  [r5(k)  +T^(k)]  p-'.-'(k) 

corr 

+  E  [r:;d(k,q)  +  rj',{k,q)]  P-’-'ik  +  q).  (8) 

q 

The  real  part  of  T"  describes  the  dephasing  of  the  polarization  due  to  scattering  of  carriers  in 
the  bands  v.  The  diagonal  contribution.  Re  Tj,  leads  to  a  momentum-dependent  polarization 
decay  rate  whereas  the  non-diagonal  contribution.  Re  mixes  the  polarizations  of  various 
k-states.  Without  higher-order  polarization  terms  [8,11],  which  can  be  neglected  in  the 
regime  where  a  weak  probe  field  is  treated,  we  find  Re  ra(k)  =  S[^(k)  +  S^j(k)  where 

=  i:r„(k)  [1  -  r(k)i  -  sj.,(k)  /"(k)  (9) 

^  corr 

determines  the  redistribution  of  the  carrier  occupation  probability  /*".  Hence  the  diagonal 


dt 


damping  rate  contains  the  sum  of  in  and  out-scattering  rates  in  the  carrier  dynamics.  For  the 
considered  quasi-equilibrium  situation  of  the  carrier  system,  the  detailed  balance  condition 
leads  to  ^/*'(k)|corr=0  whereas  in  and  out-scattering  results  in  a  large  damping  rate  Re  Fa 
which  corresponds  to  a  decay  time  on  the  order  of  100  fs.  In  Eq.  (8),  this  diagonal  damping 
rate  is  compensated  to  a  large  exten^tby  the  non-diagonal  damping  rate  Re  Fnd-  In  Ref.  [9] 
it  has  been  shown,  that  this  compensation  has  to  be  considered  in  order  to  obtain  the  correct 
carrier  generation  rate  for  intense  optical  interband  excitation.  Also  the  excitonic  damping 
and  the  lineshape  of  the  gain  [10]  can  be  properly  described  only  by  including  both 
diagonal  and  non-diagonal  damping. 

On  the  other  hand,  the  imaginary  part  of  F"  combines  with  the  mean-field  contributions, 
described  by  Eqs.  (4)  and  (5),  by  adding  screening  contributions  to  the  renormalized  carrier 
energy  and  §a^the  renormalized  Rabi  energy.  Note,  however,  that  in  general  the  correlation 
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contributions  cannot  be  written  as  corrections  to  single  particle  properties. 

In  the  following  we  consider  the  influence  of  carrier-carrier  Coulomb  scattering  as  well  as 
carrier-phonon  scattering.  In  second  Born  approximation,  where  terms  up  to  the  quadratic 
order  in  the  screened  Coulomb  potential  are  included,  the  carrier-carrier  scattering  leads  to 
the  diagonal  contribution 

rf ')(k)  =  i  E  E  {2|»',l"  -  WV+q-k')  9(4  -  ei:+,  +  4  -  4-q) 

^  k',q  I/' 

X  [4+q  (l  -  fi)  /i^-q  +  (l  -  /i(+q)  fi  (l  "  4-,)]  ,  (10) 

whereas  the  non-diagonal  scattering  rate  follows  from 

rS“’(k.  q)  =  T  E  E  w'k+,-k')  9(-4 + 4+q  -  4  +  4-q) 

^  k'  1/ 

X  [4  /i?  (1  -  +  (1  -  4)  (1  -  /£:)  /i?-,] .  (11) 

In  the  subband  matrix  element  of  the  screened  2d-Coulomb  potential,  Wq  =  Ki/cq,  the 
influence  of  excitation-induced  screening  is  described  by  the  dielectric  function  Cq  which  will 
be  calculated  using  the  static  Lindhard  formula  [1^.  Screening  contributions  of  the  crystal 
lattice  (including  phonons)  and  nonresonant  transitions  are  included  through  a  background 
dielectric  constant  €4  in  Eq.  (7).  With  =  c©  background  contributions  of  phonons  are 
taken  into  account.  However,  as  part  of  the  phonons  and  their  interaction  with  the  carriers 
will  be  treated  explicitly,  this  dielectric  constant  has  to  be  reduced  [12].  Hence,  when  we 
consider  correlation  contributions  due  to  carrier-phonon  interaction,  we  use  a  background 
dielectric  constant  c*  =  Coo  which  does  not  include  the  influence  of  phonons.  In  this  limiting 
case,  the  background  phonon  screening  is  slightly  underestimated. 

For  the  electron-phonon  scattering,  the  diagonal  rates  in  Eq.  (8)  are  given  by 

^(k)  =  —  7q  { ^{1  —  /k_q)  nr,o  +  /k_q  (1  +  n.;)o)]  ~  ^k-q  +  ^r,o) 

fi  q 

+  [(1  -  /k-q)  (1  +  n«)  +  /k-q  n/.o]  9(4  -  4-,  -  ,  (12) 

and  the  non-diagonal  rates  are 
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q)  =  T  H  7q  {[(1  -  /k)  njx)  +  /k  (1  +  nj,o)]  -  ej  +  Juvr^o) 


a=u,t/ 


+  [(1  ~  /k)  (1  +  ^TX))  +  /k  l^w]  di^i-q  -  Ck  “  f>^Lo)\  >  (13) 


where  tu&i,o  is  the  energy  of  the  LO-phonons  and  ji  is  the  Prohlich  interaction  matrix 


element  of  3d  LO-phonons  and  confined  electrons  with  in-plane  momentum  q, 


7q 


-^jdq. 


A-keo  €(|qp+<,|)  (|q|2  +  ql) 


(  _1 _ ^T.O 

\  Coo  €o  /  2 


(14) 


For  the  distribution  function  of  phonons,  nio,  we  use  a  Bose-Einstein  function  containing 
the  lattice  temperature. 

The  derivation  of  the  scattering  rates  requires  many  body  techniques  such  as  Green’s 
functions  or  projection  formalism  in  density  matrix  theory.  Here,  within  the  assumption  that 
gain  spectra  at  room  temperature  can  be  described  by  using  one-particle  correlations,  higher 
order  correlation  functions  have  been  factorized.  In  the  Green’s  functions  technique,  the 
random-phase  approximation  (RPA)  as  well  as  the  first  Coulomb  vertex  contribution  have  to 
be  considered  to  obtain  all  scattering  contributions  in  Eqs.  (10)  and  (11)  which  are  quadratic 
in  the  screened  Coulomb  potential  W.  Correlation  contributions  have  a  complicated  time 
dependence,  which  includes  memory  effects.  In  Eqs.  (10)  -  (13),  this  time-dependence  has 
been  approximated  within  an  adiabatic  treatment  that  leads  to  a  generalized  Heitler-Zeta 
function. 


9(^)  = 


e  +  i'j 


(15) 


This  Heitler-Zeta  function  includes  the  effective  quasi-particle  broadening  7  which  will  be 
treated  as  a  small  constant.  While  7  is  a  property  of  correlated  carriers,  we  do  not  include 
any  phenomenolagical  polarization  decay  time. 

At  the  end  of  this  section,  we  outline  some  properties  of  the  discussed  theory  for  a  4  nm 
CdZnSe  quantum  well  (see  next  section).  The  dependence  of  our  results  on  the  effective 
quasi-particle  broadening  7  is  shown  in  Fig.  1  where,  for  simplicity,  only  a  single  valence 
band  has  been  considered.  While  the  diagonal  scattering  rate  is  reduced  with  increasing 
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7,  the  compensation  of  diagonal  and  non-diagonal  rates  results  in  a  spectrum  which  is 
basically  independent  on  7.  In  Fig.  2  we  study  the  influence  of  various  approximations 
for  the  polarization  dephasing  on  the  absorption  spectrum.  In  the  simplest  approximation, 
correlation  contributions  in  Eq.  (2)  are  included  phenomenologically.  Polarization  dephasing, 
which  is  described  by  the  real  part  of  Eq.  (8),  has  often  been  approximated  by  a  constant  T2 
time.  Screening  corrections  to  the  Hartree-Fock  terms,  Eqs.  (4)  and  (5),  which  follow  from 
the  imaginary  part  of  Eq.  (8),  have  been  approximated  by  replacing  V  hy  W  and  adding 
the  Coulomb-hole  term  in  Eq.  (4)  [12].  The  resulting  absorption  spectra  are  shown  in 
Fig.  2a.  With  increasing  carrier  density,  we  obtain  a  strong  artificial  red  shift  of  the  exciton 
resonance  (at  moderate  densities)  and  of  the  gain  peak  (at  high  densities).  As  a  next  step, 
we  consider  correlation  contributions  due  to  carrier-carrier  scattering.  In  Fig.  2b,  only  the 
complex  diagonal  rate  rj^“^(k)  has  been  used  together  with  Eqs.  (4)-(8),  whereas  in  Fig.  2c 
also  the  non-diagonal  rate  rnd“^(k,q)  has  been  included.  When  only  the  diagonal  rate  is 
considered,  dephasing  is  clearly  overestimated  which  can  be  seen  from  the  strong  broadening 
of  the  exciton  resonances.  Only  if  both  diagonal  and  non-diagonal  rates  are  included,  the 
exciton  broadening  is  reduced  to  reasonable  values  and,  in  agreement  with  experimental 
observations,  almost  no  shift  of  the  exciton  resonance  is  obtained  with  increasing  carrier 
density.  In  Figs.  2d  and  e,  we  directly  compare  the  approximations  of  Figs.  2a-c  for  carrier 
densities  leading  to  gain.  In  the  phenomenological  model  (dotted  line)  as  well  as  in  the  pure 
dephasing  limit  (dashed  line),  the  bandgap  renormalization  is  drastically  overestimated, 
the  lineshape  of  the  gain  is  modified  and  unphysical  absorption  below  the  renormalized 
bandedge  is  obtained  in  comparison  to  the  full  model  (solid  line).  Especially  the  absorption 
energetically  below  the  gain  region,  which  occurs  in  simplified  gain  calculations,  is  a  signature 
of  overestimated  dephasing  and  incorrect  lineshape.  The  resulting  strong  broadening  of 
carrier  lineshape  functions  leads  to  an  unphysical  admixture  of  non-inverted  states  high 
above  the  band  edge. 

In  conclusion,  the  large  red  shift  of  the  exciton  at  low  carrier  densities,  the  overestimated 
band-gap  shift  at  higher  densities  as  well  as  the  large  damping  are  artifacts  of  simpler 
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models  which  are  in  clear  contradiction  to  experimental  results.  Hence,  a  pure  dephasing 
approximation  is  not  appropriate  and  only  the  full  model  can  reproduce  the  experimental 
results. 

For  a  realistic  description  of  quantum  wells  on  the  basis  of  II- VI  compounds,  additional 
excitation  independent  inhomogeneous  broadening  has  to  be  considered  in  addition  to  the 
excitation  dependent  carrier-carrier  and  carrier-phonon  interaction.  In  particular,  small 
spatial  variations  of  the  concentration  of  the  quantum-well  and  barrier  materials  as  well  as 
roughness  of  the  well/barrier  interfaces  give  rise  to  inhoinogenous  broadening.  Correspond¬ 
ingly  the  calculated  results  are  more  realistic  when  the  homogeneously  broadened  spectrum 
ochom{<^)  is  convoluted  with  a  Gaussian  distribution  Q{u})  of  given  spectral  width  according 
to 

(Xinhi^o)  =  f  dco'  Oihom{<^  “  w')  (16) 

Figure  3  shows  the  influence  of  increasing  additional  inhomogenous  broadening  on  the  TE- 
spectrum  of  a  4  nm  CdZnSe  quantum  well  (see  next  section).  Note  that  the  influence 
of  inhomogeneous  broadening  is  stronger  for  smaller  carrier  densities  where  the  spectrum 
contains  sharp  features.  With  increasing  broadening,  the  lineshape  is  modified  and  the  gain 
maxima  can  shift  several  meV.  In  the  following  results,  we  have  included  an  inhomogenous 
broadening  of  10  meV. 

III.  SPECTRA  OF  II- VI  QUANTUM  WELLS 

The  microscopic  model  discussed  in  the  previous  section  will  be  used  to  compute  the 
density-dependent  absorption  and  gain  spectra  for  CdZnSe  quantum  wells  between  15nm 
ZnSe  barriers  and  ZnSSe  cladding  layers  (25%  Cd  and  6%  S).  Two  samples  with  a  quantum- 
well  width  w=4  nm  (QWl)  and  w=7  nm  (QW2)  are  compared. 

Using  the  envelope  approximation  method  the  energy  gaps  and  the  strain  induced 
shifts  for  the  heavy  and  light  hole  states  are  calculated.  The  conduction  bands  are  always 


9 


assumed  to  be  parabolic.  Neglecting  the  split  off  band,  the  hole  states  are  coupled  via  a 
4x4  Luttinger  Hamiltonian.  The  diagonalization  of  this  system  leads  to  the  valence  bands 
eu>0^)  and  the  dipol  matrix  elements  All  those  bands  are  considered  where  the  carriers 
are  confined  in  the  quantum  well  for  at  least  small  wavevectors  k.  For  the  4  nm  and  7 
nm  quantum  wells,  the  calculated  band  structure  is  shown  in  Fig.  4.  In  the  4  nm  (7  nm) 
quantum  well,  one  (two)  conductions  band  and  three  (five)  valence  bands  are  confined.  The 
compressive  strain  of  the  CdZnSe  wells  between  the  ZnSe  barriers  produces  a  heavy-hole 
light-hole  splitting.  For  the  7  nm  quantum  well,  the  top  three  valence  bands  are  heavy-hole 
like  close  to  the  band  edge  so  that  their  coupling  to  the  TM-mode  is  weak.  Hence  the 
exdton  and  gain  corresponding  to  these  transitions  is  small.  The  main  contribution  to  the 
TM-mode  follows  from  tansitions  between  the  fourth  and  fifth  valence  band  and  the  first 
conduction  band. 

The  TE-  and  TM-spectra  of  the  4  nm  and  7  nm  quantum  wells  are  plotted  in  Figs.  5 
and  6,  respectively.  The  shown  carrier  densities  cover  the  region  from  no  inversion  between 
any  bands  to  inversion  between  most  of  the  confined  bands  for  small  wavevectors  k  [15]. 

The  low-density  TE-spectra  of  the  4  nm  quantum  well  (left  part  in  Fig.  5)  show  two  strong 
excitonic  absorption  lines  corresponding  to  the  1-1  and  2-2  transition.  Here  transitions  are 
labeled  i-j  with  i  (j)  refering  to  the  conduction  (valence)  band  involved.  For  our  system,  the 
second  TE-absorption  (2-2)  line  appears  at  the  same  spectral  position  as  the  TM-absorption 
(1-4)  line;  the  differences  of  the  band  gaps  for  these  transitions  is  less  than  2meV.  These 
excitonic  absorption  lines  stay  at  their  spectral  position  when  the  carrier  density  increases; 
after  exciton  bleaching  the  transition  develops  gain.  For  the  highest  carrier  density,  the  TE- 
spectra  of  the  4  nm  quantum  well  exhibit  a  shoulder  due  to  inversion  of  the  2-2  transition.  In 
addition,  the  left  part  of  Fig.  5  shows  that  the  absorption  above  the  renormalized  band  edge 
does  not  decrease  when  the  carrier  density  increases;  the  high-energy  part  of  the  spectra  for 
large  carrier  densities  is  only  shifted  as  a  consequence  of  bandgap  renormalization  due  to 
Coulomb  interaction. 

If  the  quantum  well  is  wider,  more  bands  are  confined.  This  is  illustrated  in  the  right 


part  of  Figs.  5  and  6.  Due  to  the  second  conduction  band  in  QW2,  we  obtain  an  additional 
excitonic  transition  in  the  TE-spectra  while  the  TM-absorption  remains  similar  to  that  of 
QWl.  As  the  carrier  density  is  distributed  over  more  bands  in  QW2,  the  inversion  becomes 
smaller  and  the  gain  is  reduced.  In  addition,  the  polarization  decay  is  increased  for  QW2 
at  low  plasma  densities  which  leads  to  a  broader  excitonic  absorption  line. 

The  band  offset,  i.e.,  the  splitting  of  the  confinement  energy  between  the  conduction  and 
the  valence  band,  cannot  be  determined  unambiguously  from  experiments.  In  the  previous 
calculations,  a  band  offset  of  60%  for  the  conduction  band  has  been  assumed.  The  influence 
of  band  offset  changes  on  our  results  is  studied  in  Figs.  7  and  8  for  the  TE-  and  TM-spectra, 

respectively.  For  a  conduction  band  offset  between  30%  and  60%,  a  single  sub-band  appears 

- - - 

in  the  finite-hight  quantum  well  confinement  potential  of  the  electrons  using  the  above 
discussed  material  composition  and  4  nm  well  width  (QWl).  However,  the  confinement 
potential  of  holes  leads  to  four  sub-bands  within  the  quantum  well  for  30%  and  40%  offset 
and  three  sub-bands  when  50%  and  60%  offset  are  considered.  The  TE-spectra  for  various 
band-offset  values  are  shown  in  Fig.  7.  The  lineshape  of  the  spectra  at  high  densities  is 
almost  independent  of  the  band  offset;  only  the  magnitude  of  the  gain  is  slightly  changed. 
Also  the  low-density  spectra  are  similar,  only  the  weak  excitonic  resonance  of  the  fourth 
hole  subband  is  missing  for  a  band  offset  >  50%.  With  increasing  band  offset,  the  envelope 
function  of  the  light  hole  becomes  less  confined  in  the  quantum  well  which  leads  to  an 
reduced  effective  TM-dipol'^oupling  with  the  conduction  band.  For  that  reason  both  the 
TM-absorption  for  low  plasma  densities  and  TM-gain  at  higher  densities  are  reduced  with 
increasing  band  offset  as  shown  in  Fig.  8. 

As  another  application  of  our  theory,  we  study  the  influence  of  compressive  versus  tensile 
strain  on  the  quantum- well  spectra  of  II- VI  compounds.  We  investigate  a  4  nm  Zno.8Cdo.2Se 
quantum  well  with  a  lattice  constant  of  5.75  A  between  ZnMgSSe  barriers  having  a  lattice 
constant  of  either  5.65  A  (QW3)  or  5.83  A  (QW4)  [16].  The  corresponding  band  structure 
is  shown  in  Fig.  9.  In  QW3,  the  first  three  valence  bands  are  heavy-hole  like  at  the  zone 
center  while  in  QW4  the  tensile  strain  causes  a  top  valence  band  with  light-hole  character. 
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Therefore,  the  TE-gain  strongly  exceeds  the  TM-gain  for  compressive  strain  whereas  for 
tensile  strain  the  TM-gain  is  dominant  (Fig.  10).  The  carrier  density  dependence  of  the 
TE-spectrum  for  QW3  is  shown  in  Fig.  11.  At  intermediate  densities,  a  situation  can  be 
realized  where  a  small  gain  exists  while  the  excitonic  enhancement  is  still  present.  A  similar 
situation  has  recently  been  observed  in  ZnCdSe/ZnSSe/ZnMgSSe  quantum  wells  [17]. 


IV.  CONCLUSION 


In  summary,  the  gain  spectra  of  II- VI  quantum- well  materials  have  been  investigated 
within  a  microscopic  plasma  theory  which  is  based  on  kinetic  equations  for  an  interacting 
electron-hole  system  in  a  multi-band  semiconductor.  These  equations  include  many-body 
effects  such  as  intraband  as  well  as  interband  Coulomb  correlations  leading  to  carrier  scatter¬ 
ing  and  excitonic  transitions.  In  addition,  the  influence  of  LO-phonons  has  been  considered. 
The  single-particle  energies  and  the  dipole  coupling  have  computed  using  Luttinger-Kohn 
theory  in  envelope  function  approximation. 

The  model  has  been  used  for  calculating  the  excitation  density  dependent  absorption 
and  gain  spectra  of  several  different  examples.  The  possibility  to  optimize  the  gain  spectra 
in  terms  of  well-width  and  composition  has  been  studied  for  II- VI  compounds. 
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FIGURES 


FIG.  1.  Complex  scattering  rate  for  electrons  (e)  and  holes  (h)  and  absorption  spec¬ 

trum  for  a  two-band  model  with  carrier  density  2  x  10^^  cm~^  and  various  quasi-particle  broaden¬ 
ing  7. 

FIG.  2.  Influence  of  the  polarization-dephasing  model  on  the  linear  absorption  spectrum  of 
a  4  nm  CdZnSe  quantum  well.  Within  a  phenomenological  model  (a),  a  constant  polarization 
decay  rate  T2=100  fs  and  screened  Hartree-Fock  contributions  have  been  considered.  Using  a 
microscopic  model  for  correlation  contributions  due  to  carrier-carrier  scattering,  we  compare  the 
diagonal  dephasing  (b)  with  diagonal  and  non-diagonal  dephasing  (c).  The  carrier  densities  are 
0.1,  0.5,  1,  2,  4  xlO^^  cm~^  (from  top  to  bottom).  In  (d)  and  (e)  we  directly  compare  the 
phenomenological  model  (dotted  line)  with  the  diagonal  dephasing  (dashed  line)  and  diagonal  + 
non-diagonal  dephasing  (solid  line)  for  a  fixed  carrier  density. 

FIG.  3.  Absorption  spectrum  of  a  4  nm  quantum  well  for  a  carrier  density  2  x  10^^  cm“^  and 
various  inhomogeneous  broadening. 

FIG.  4.  Energy  bands  for  a  4  nm  and  7  nm  CdZnSe/ZnSe  quantum  well.  The  offset  energy  Ei 
is  2.52  eV. 

FIG.  5.  Comparison  of  the  TE-spectra  of  a  CdZnSe/ZnSe  quantum  well  with  4  nm  well  width 
(QWl)  and  7  nm  well  width  (QW2).  The  carrier  densities  are  0.5,  1,  2,  3,  4,  5,  6,  7,  8  x  10^^  cm“^ 
(from  top  to  bottom).  An  additional  inhomogenous  broadening  of  10  meV  is  taken  into  account. 
The  bottom  figures  show  the  same  on  a  larger  scale. 

FIG.  6.  Same  as  Fig.  5  for  the  TM-spectra. 

FIG.  7.  TEl-spectra  for  QWl  where  different  band  offsets  are  assumed.  The  densities  are  similar 
as  in  Fig.(6).  For  30%  or  40%  conduction-band  offset,  four  valence  bcinds  are  confined  in  the  well 
and  three  valence  bands  for  higher  offsets. 
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FIG.  8.  Same  as  Fig.  7  for  the  TM-spectra.  By  increasing  the  conduction-band  offset,  i.e., 
reducing  the  potential  well  for  the  holes,  the  light  hole  becomes  less  confined  and  therefore  the 
TM-coupling  of  the  conduction  bands  with  the  confined  valence  bands  decreases. 

FIG.  9.  Energy  bands  for  a  4  nm  CdZnSe/ZnMgSSe  quantum  well  with  compressive  (QW3) 
and  tensile  (QW4)  strain.  The  offset  energy  E2  is  2.68  eV. 

FIG.  10.  TE-  and  TM-spectra  for  a  4  nm  CdZnSe/ZnMgSSe  quantum  well  with  compressive 
(QW3)  and  tensile  (QW4)  strain  at  carrier  density  6  xlO^^  cm“^.  The  bottom  figures  show  the 
seime  on  a  larger  scale. 

FIG.  11.  TE-spectra  for  a  4  nm  CdZnSe/ZnMgSSe  quantum  well  with  compressive  strain 
{QW3)  for  carrier  densities  1.5,  1.8,  2.0,  2.2,  2.5  and  4.0  xlO^^  cm“^  (from  top  to  bottom). 
The  bottom  figure  shows  the  saune  on  a  larger  scale. 
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